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We study molecular hydrogen formation in and on solids. We construct a model with surface sites 
and bulk sites capable of describing (1) the motion and exchange of H and H 2 between surface and 
bulk, (2) the recombination of H and dissociation of H 2 in and on the solid, and (3) the injection of 
H from the gas phase and the loss of H and H 2 from the solid. The basic physical processes include 
thermally activated reactions, collisionally induced reactions, and tunneling reactions. 

Our main application is to the astronomical problem of H 2 formation on grains in space but the 
model has more general applicability. We investigate the steady-state H and H 2 concentrations in 
and on the solid when gas phase atoms or ions stick to the surface or penetrate the body of the grain. 
The model identifies ranges of physical parameters for which the solid becomes saturated (surface 
and/or bulk) with particles (H and/or H 2 ) and facilitates the calculation of the efficiency of molecule 
formation (the fraction of the gas phase atoms that leave as molecules). These solutions are highly 
degenerate in the sense that they depend only on a small number of dimensionless parameters. We 
find that a variety of recombination pathways operate under a broad range of conditions. 

As an example we study H 2 formation in and on carbon grains. If molecules form from H on 
the surface by quantum tunneling alone, efficient transformation of incoming atoms takes place at 
grain temperatures less than 75 — 100 K for 1 < (n(H)/cm _3 )(e/eV) 1 ^ 2 < 10' where n(H) and 
e are the number density and energy of incident particles. When incident particles are energetic 
enough to penetrate the surface, the bulk of the grain will be saturated by H and/or H 2 under most 
circumstances. Additional molecule formation pathways (recombination on the surface or in the 
bulk) increase the efficiency of the transformation but only rapid bulk to surface exchange of H 2 
can alter the conclusion regarding saturation. Saturation can lead to fundamental changes in the 
parameters that describe the grain. 
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I. INTRODUCTION 


The interaction of hydrogen with solids is a topic of considerable interest in astrophysics, as well as other applied 
fields (materials research in fusion devices, thin film growth and surface etching of semiconductors, and adsorption 
and desorption reaction dynamics in catalysis). There is likewise a great deal of interest in exploring the fundamental 
physical processes that govern the interaction between gas atoms/ions and the solid surface/bulk. In this paper we 
develop a model to describe the behavior of H and H 2 in and on solids. The model is phenomenological - it identifies 
and parameterizes a set of reaction pathways thought to play important roles in determining the concentrations of H 
and H 2 in and on the solid. 

Our main application is the study of efficiency of molecule formation on and release from astronomical grains. 
The possibility that grain surfaces catalyze the formation of H 2 inside cold atomic/molecular clouds has long been 
recognizedB. Under quiescent conditions, the chemisorption sites are filled and hydrogen is weakly bound by physical 
adsorption. In this scenario, low grain temperatures are necessary for the catalysis else the H atoms evaporate before 
they can recombine. The gas temperature must also be low else the neutral gas atoms rebound from the grain 
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surface. Observations show that molecular hydrogen emission is excited d 
studiedj'egions include the Orion Molecular Cloud (Orion OMC-1, e.g. refs.l , 

CygnuaS), various Herbig-Haro objects (HH 7-11, e.g., refs.Q’ta, Seyfert galactic nuclei (NGC1275H) and interacting 
starburst galaxies (NGC3690-IC694u). In each cited example the emission has been linked to shock heating, however, 
shock models often run into a common difficulty, finding a way to account for the emitting molecular hydrogen under 
decidedly non-quiescent conditions. 

The main motivation is to explore new and different pathways for H 2 formation that may operate under more 
extreme conditions than hitherto considered, e.g. in the vicinity of a 100 km s _1 shock. Several possibilities are 
of interest: (1) Energetic protons can penetrate the grain lattice and diffuse within. H can recombine within the 
grain or diffuse to the surface and recombine. (2) Energetic particles can drive endothermic molecule formation 
reactions and also eject and/or dissociate H containing species from the grain. (3) Energetic particles can sputter 
the surface and generate a large surface density of pits and edges, sites which can act as points of increased physical 
adsorption. We explore the important physical parameters that regulate molecule formation in the first two scenarios. 
These possibilities may turn out to be of astrophysical importance if they can provide means for rapidly reforming a 
small amount of molecular hydrogen in gas that has recently been shocked. If fast shocks radiate significantly more 
molecular line radiation than is currently thought, re-interpretation of a large body of work may be necessary. 

Motivated by these considerations, we explore grain mediated molecule formation under a wide range of conditions 
including those that would be found in the vicinity of a shock. We describe the most important physical parameters, 
we review which parameters have been determined by laboratory experiments and astronomical observations, and 
we construct models that illustrate the different qualitative outcomes for molecule formation. Although our analysis 
is motivated primarily by astrophysical considerations, it is sufficiently general to be applied to other systems. For 
example, modeling high fluence hydrogen implantation in solids may follow the general strategy presented in this 
paper. 

In §jnj we provide a brief overview of the typical energy barriers of interest in the formation pathways and, for 
comparison, the particle energy scales expected behind an astrophysical shock. We describe the typical grain tem¬ 
perature and note the wide range of grain composition and size with which we are concerned. In § III we outline the 
molecule formation pathways and identify key parameters. We briefly review the following topics: the transmission 
and reflection of incoming gas particles, the stopping ranges, the possible ordering of energy levels between the bulk, 
surface and vacuum, the H-H pairwise interaction in and on the solid, and the diffusion of H in and on the solid. 
Where necessary, we adopt a more or less phenomenological description of rate coefficients for the succeeding analysis. 
In ' JFv| we formulate the model, including the basic equations and the conservation laws. We focus on a grain with one 
type of surface site and one type of bulk site. H 2 formation is calculated in three successively more general models: 
(i) pairwise recombination involving H 2 ejection from the surface, (ii) plus recombination and binding of H 2 on the 
surface, (iii) plus recombination and binding in the bulk. In §|v| we apply our results to graphite grains. In §VI we 
summarize. 


II. CHARACTERISTIC POSTSHOCK ENERGY SCALES 

Three key endothermic processes - ion implantation in the solid, sputtering of the surface atoms and chemical 
reactions driven by energetic particles - occur in the molecule formation pathways. As we discuss in more detail 
below, the threshold for implantation is expected to be ~ 10 eV, but has not been well-characterized for astrophysical 
grain material. The sputtering thresholds for H (He) range from ~ 25 — 60 eV (~ 10 — 16 eV) for refractory grain 
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material (graphite, silicate and irorJl3). The strength of a typical H bond is 1 — 5 eV (chemical) or ~ 0.05 — 0.5 eV 
(physical absorption, El). The equipartition thermal energy of particles behind a 100 km s -1 shock is ~ 10 eV and 
several factors may increase the energy of gas-grain-cpllisions. As the grains of size a interact with the plasma they 
typically acquire a negative charge Q ~ — 2.5a/cT/eE3 (a number of important additional effects are noted in ref.Ej) 
so that protons strike grains with a total energy of ~ 35 eV. In addition, grains of sufficiently large size (a > ICC. 5 
cm) are betatron accelerated behind a strong radiative shock to velocities ~ 3v s , where v s is the shock velocitycl. 
From the shock front to the recombination region, the grains interact with H + , He + and He ++ with typical energies 
from 1-300 eV. The characteristic timescale for the gas to cool to 10 4 K is t coo i ~ 2 x 10 10 i; 5 7 2 /no sEJ, where no is the 
preshock H-nuclei density and v s t = v s /10 7 cm s — . 

Observations show that interstellar grains span a range in size at least as wide as 5 x 10 -7 < a < 2 x 10 -5 cm 
and have absorption/emission features-typical of silicate (Si-0 stretching and O-Si-O bending modes), graphite and 
hydrocarbons (C-H stretching modes)e3. 

The ratio of the grain temperature to the hydrogen binding energy is an essential parameter in all molecule 
formation schemes; it regulates the rate of thermally activated diffusion, recombination and evaporation of atoms 
and molecules from the surface. The grain temperature is largely governed by the balance of heating of Lyman- 
alpha radiation trapped in the vicinity of the shock front with cooling by infrared emission. Silicate grains reach 
T gr « 320K [(?r o /10 6 )uf 7 /a_5]°' 2 and the results for graphite grains are comparable] 44 !. Note that at low densities, 
T gr is much less, e.g. T gr ss 20K °’ 2 and quantum mechanical tunneling competes with the mechanism of 

thermal activation in the processes listed above. 


III. KEY PARAMETERS 


A. Outline for Molecule Formation 


Here we provide a schematic overview of the formation mechanism: 

1. H or H + impinges on a grain; if H + , it is neutralized by electron transfer as it approaches. 

2. A collision takes place and the atom rebounds, sticks or penetrates the surface. The lattice may be left intact, 
sputtered or damaged (i.e. defects introduced). If the grain surface includes H physically or chemically bound (X-H), 
then the incoming particle may drive a “pickup” reaction of the form X-H + H —> X + H 2 . 

3. Within the grain, the atom diffuses from site to site. It seeks out traps (vacancies, interstitials, etc.) where its 
binding energy is greatest. An equilibrium distribution is eventually formed in which the tightest binding sites are 
preferentially occupied. 

4. When a new atom enters a grain it explores the grain interior by diffusion. Either it reaches the surface or it 
finds another atom, overcomes the activation barrier and forms H 2 . 

5. Bulk H 2 diffuses out of the grain, or remains in situ until it can escape directly. This is possible when damage 
of the lattice (by implantation, sputtering, radiation damage or grain-grain collisions) has accumulated to the point 
that passageways to the surface form. 

6. If H 2 formation in the bulk is energetically unfavorable or inefficient, then the newly trapped atom diffuses to 
the surface where it recombines with another atom to form an H 2 molecule which immediately escapes from the grain. 
If the grain is too cold for-surface recombination to proceed by thermal activation then quantum tunneling may allow 
recombination to proccecQ 

7. In the case in which the-ambient atoms are not energetic enough to penetrate the grain, the H 2 formation 
scenario suitable for cold cloudsEl takes place, but at possibly higher grain temperature and in the presence of enhanced 
binding sites (referred to as “semi chemisorption” in that paper). Atoms fill up the tightest binding sites, until the 
recombination rate plus the atom’s evaporation rate balances the flux of incoming particles. 

The chemical network is governed by a number of key parameters: the mean interception time of H nuclei by the 
grain, the diffusion time of H and H 2 inside the grain, the surface and bulk binding energies for H and H 2 and the 
energy barriers for H 2 formation. In addition to normal sites in and on the grain, there are likely to be impurity sites 
whose number and characteristic binding energy for H and H 2 are important. 

In the succeeding subsections we review some of the information relevant to the determination of these parameters 
and, ultimately, to our characterization of the rate coefficients for the model. In section IV, we examine in detail how 
the efficiency of molecular hydrogen formation depends on such quantities. 
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B. Transmission and Reflection of Low Energy Particles 


At low energies (E <100 eV) the ab initio calculation of the transmission and reflection coefficients of particles from 
a target surface is a formidable task. The results are^otimportance in plasma wall experiments, and such calculations 
have been the subject of many recent investigationsiJ’H. Many-body effects play an especially important role. Firstly, 
in contrast to the scattering of high-energy particles, which can be well approximated as a pure binary collision between 
the projectile and a surface atom, the scattering of low-energy particles involves simultaneous interaction with several 
nearby atoms, owing to the relatively long-range force between the ion and each surface atom (e.g., refs.tu). As a result, 
the scattering of low energy ions is sensitive to the details of the ion-surface interaction potential. Secondly, in the 
case of scatteringJi'pm metal surfaces, collective effects produce a significant attractive force, commonly represented 
as an image forceE3c2l, which plays a major role in determining the transmission coefficient at low energies. 

Ignoring chemical reactions, there are three possible outcomes for a fast particle incident on a surface: the incident 
particle may stick to the surface of the target material, penetrate below the surface or reflect back into space. The. 
typical energy at which each outcome is maximized varies with the-incident ion and the target material. Maselcfl 
reviews the status of trapping and sticking on solid surfaces. BaskesLj performed simulations of hydrogen reflection 
from a clean nickel surface, using the Embedded Atom Method to handle the many-body interactions in a self- 
consistent manner. For normal incidence, he found that the reflection coefficient peaks at an energy of order 6 eV 
where about 90% of the incident ions are reflected. Well above 6 eV the ions are energetic enough to penetrate the 
surface, and below it the ions are trapped, owing to the attractive image force which tends to increase the energy loss 
of the scattered ion tp-the surface. Similar theoretical results were obtained by Eckstein & Biersacklifl who employed a 
modified TRIM codecil to calculate hydrogen reflection. They treated the surface binding energy as a free parameter 
and showed that at energies well above a few eV the reflection coefficient is independent of the surface potential, but 
at < a few eV, depending on the surface binding energy, becomes sensitive to it. In experiments similar behavior has 
been observed for scattering of Na + from CuEi 

Less is known about the reflection and sticking of H to non-metals. The sticking probability of atomic-hydrogen 
on graphite_at low (sub-eV) energies has been inferred from measurements to be at least a few percento TRIM 
simulationain of reflection of low energy H atoms from carbon give a maximum reflection coefficient smaller than from 
Ni and W, and a slightly higher energy at the maximum. These differences are mainly due to the difference in target 
mass. For a surface binding energy of 1 eV about 50% of the H atoms are found to be reflected from the carbon 
surface at the peak energy of 4 eV. 

Let T(E , a) be the fraction of incident particles that stick to or penetrate and do not exit a grain of size a at a given 
energy E. (Hence T ~ 0 for particles whose energy-dependent stopping range exceeds the grain size; equivalently, 
T(E,a) ~ 0 for E > E max (a) when E max is the maximum energy particle stopped by the grain.) Let a g be the 
grain’s geometric cross section for intercepting particles (including the focusing effect of Coulomb scattering), let v g 
be the grain’s velocity and let f(v) be the particle distribution function. Then the mean rate for particle interception 
by the grain is 

tin = J d 3 vf(v)a g \v-v g \T(E,a). (1) 


C. Stopping Ranges of Low Energy Ions in Solids 

Energetic particles that penetrate into the grain will experience energy loss due to nuclear scattering and electronic 
stopping, i.e., energy transfer from the ion to the target nuclei and electrons (for a detailed account of the theory 
of ion stopping in solids see ref.Ed) . In order to be trapped inside the grain, the particle’s stopping range must not 
exceed thu characteristic grain size. At low energies, the stopping power is dominated by nuclear scatteringpl (see 
also refs £3 for reviews on the stopping power of an electron gas). 

The energy transferred in a pure binary collision from an incident projectile having energy Eo and mass to, to a 
target particle of mass M, is 


aH, = 4 ^(TfW si “1' 


( 2 ) 


where /i = m/M is the mass ratio, and 9 is the center-of-mass scattering angle. The lab frame scattering angle, ip, 
satisfies 
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tan if) = 


sin 9 

cos 9 + n' 


( 3 ) 


For a proton scattering off heavy target nuclei (/j, « 1), Eq. || shows that the fraction of energy lost by the intruder 
in one collision is at most of order /r. Consequently, the number of collisions required to slow the incident proton from 
Eq to final energy Ef will be ~ ln(Eo/ Ef) //z, and the corresponding stopping range, assuming that the ion undergoes 
a random walk with fixed step size, varies like /z -1 / 2 . However, as already mentioned in § III B, at low energies, the 


scattering of an ion from an array of atoms is considerably more complicated than pure binary collision because the 
ion interacts simultaneously with several nearby atoms. This tends to increase the effective mass of the target object, 
thereby reducing the stopping power. Moreover, screening plays an important role at these energies and needsJto 
be taken into account properly. Quantitative determination of range distributions requires numerical simulations^. 
For our purposes, however, a rough upper limit on the average stopping range should suffice. For incident energies 
Eq < 100 eV, final energy Ef ~ 1 eV, fj, ~ 1/12 (e.g. graphite grain), we estimate < 70 collisions. If the mean free 
path is not greater than the lattice spacing then the stopping range of H anil not exceed ~ 8 atomic layers. This 
upper limit is in agreement with recent experimental results reported in refHj. 

We assume the grains are large enough to stop the particles that penetrate and we treat the system as a semi-infinite 
sample. 


D. Energy Levels 


The ground state energy of H and H 2 in the bulk, on the surface and in the vacuum are the fundamental parameters 
that govern the behavior of hydrogen with respect to solids. We denote the ground state energy of species X in the 
bulk, on the surface and in the vacuum as E [X] 1 where i = b, s, and v respectively. For future use we define 

• Chemisorption energy (for H) E c = E[H] S — ^E[H 2 \ V . 

• Solution energy (for H) E s = E[H] fa - |E[H 2 ] w . 

• Dissociation energy of H 2 (per nuclei) in bulk, on surface and in vacuum E l d = E[H] 1 — iE[H 2 ]*. 

• Embedding energy (for H) E e = Epf]" — E[H] S . 


The most basic properties of the solid are governed by the signs and sizes of these energies. We adopt as the zero 
of the energy scale the level associated with H 2 in the vacuum, i.e. E[H 2 ] 1i = 0, so it follows E^ = 2.24 eV. For an 
arbitrary grain there remain 4 quantities to specify. A simple counting shows that there are 3x4x5x6 = 360 
possible orderings of the remaining levels, so it is impractical to consider separately all relative arrangements. Even 
for H levels alone, there exist 3 x 4 = 12 possible arrangements (relative orderings) of the levels in the bulk and 
on the surface. However, general arguments regarding H 2 formation and release highlight a much smaller group of 
relevant, distinct orderings. We discuss some of these below. Note that our method of computation (section IV) 
does not depend upon an assumed ordering. Further, it should be clear that in counting the distinct possibilities 
above, we have ignored many additional energy scales that can play important roles, such as the height of barriers for 
site-to-site migration within the bulk and the height of barriers for bulk-to-surface and surface-to-vacuum transitions 
in the potential energy curves. 


1. H in bulk and surface; no H 2 

We begin by reviewing the chemisorption and solution energy scales for H. These are of considerable interest in a 
number of areas, especially fusion research on plasma-surface interactions. New experimental and theoretical . tachp 
niques have led to progress in understanding of the physics and quantitative calculation of the surface interactioriMcd. 
In many cases, the electronic wave function of the H atoms tends to form a chemical bond with the surface atoms, 
leading to a strong attractive force between the H atom and the solid surface as the atom approachs. For metals, 
the chemisorption sites are typically deeper than bulk sites (with the possible exception of deep trapping sites due to 
vacancies or defects). Eor semiconductors, bonds of a more local nature are formed and, for insulators, the situation 
is not well-understoodEd. In these cases and for graphite, the energy of chemisorption sites can lie above the energy 
of bulk sites. 

Let us consider the possibilities when the H levels satisfy the inequality {E C ,E S } < E[H]' U , which is typically the 
case for solids. The ordering implies that the surface and bulk ground states of H are energetically stable with respect 
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to H in the vacuum. The inequality drives the accumulation of H in or on the grain and the resultant high densities 
may promote H 2 formation by recombination from surface sites. Of course, H may escape by thermal evaporation, 
by collisional ejection or by molecule formation and subsequent ejection. If H and H 2 loss processes are slow then 
saturation occurs which will generally lead to changes in the solid’s properties including energy levels. Given our 
assumptions above, there exist a total of 6 distinct level orderings. 

For metals, it is typically the case that E c < and the three possible ID energy diagrams, shown in figure 0, 
have (i) 0 < E c < E s , (ii) E c < 0 < E s and (iii) E c < E s < 0. The solid line represents the schematic interaction 
potential of a single H atom with the solid. The figures illustrate the chemisorption energy E c , solution energy E s , 
and the vacuum dissociation energy of H 2 (per H) E%. The two dashed lines represent half the energy of an H 2 pair 
at a given distance from the surface. One dashed line applies to a pair at fixed interatomic separation (equal to the 
bond length of a free H 2 molecule ~ 0.74 A) and the other to the minimum energy configuration. The minimum 
energy configuration terminates near the surface where the pair separation becomes large. In general, the sign of the 
energy is defined by the direction of the corresponding arrow; positive if the arrow points upwards and negative if it 
points downwards. _. 

It is observed that (e.g,E3, and references therein) as an energetic incoming H 2 molecule approaches the solid, it 
may dissociate into H atoms that are bound to the surface (with activation energy barrier per H of E a ); in the 
inverse process, pairs of H atoms on the surface may tunnel through or thermally ascend the energy barrier and 
recombine as H 2 (with recombination energy barrier per H of E r ). In figure |l| vertical bars without arrows represent 
manifestly positive quantities. We note that quantum pair recombination from the surface is allowed only in case 
(i). Surface recombination requires thermal activation in_cases (ii-iii).[] In the case of metals, the embedding energy 
E e = £J[H]” — E c typically lies in the range 2.4 to 2.8 eVE], implying E c < 0 [case(ii) or (iii)]. _ 

The situation for carbon and other non-metals is not as well characterized. Schermann et al.Eil have reported the 
detection of H 2 formation by recombination of H atoms adsorbed on a carbon surface at temperatures between 90 and 
300 K. The detected molecules appear to be in highly excited vibrational states. At zero temperature the measured 
vibrational states of the newly formed molecules constrain the embedding energy E e < 1 eV for the H binding sites on 
the surface; consequently E c = JT[H]' U — E e > 1.24 eV. The system is so cold that thermal effects cannot significantly 
change these estimates although there remains the possibility that H may bind to more than, one type of surface site 
(for example, one physical and the other chemical). The solubility data implies E s < 0 (seeEl and references therein) 
so that E s < 0 < E c , illustrated as case (iv) in figure [j|. It is plausible that recombination proceeds via quantum 
tunneling in this experiment. The possibility of quantum pair recombination may have important implications for H 2 
formation and release which are considered below. 

Two other orderings (E s < E c < 0 and 0 < E s < E c ) are possible but we are unaware of representative materials 
for which these cases would be relevant. 


2. H and H 2 in bulk and surface 

We now review what is known about the H pair interaction on the solid surface and in the bulk. It could be 
dramatically different than that in vacuum as it depends critically on many body interactions. For most metalspas 
well as carbon, solubility data for hydrogen at sufficiently low bulk concentrations appear to follow Sieverts’ lawEn; 
that is, the bulk concentration of hydrogen is proportional to the square root of the equilibrium pressure. This 
implies that H 2 ismpt the preponderant state inside the solid, rather H (or one of its charged states) is. Some recent 
theoretical studied^ have shown that the interaction between hydrogen atoms on metal surfaces and between atoms 
in surface and subsurface sites is repulsive and short-range with interaction energies of order 0.1 - 0.4 eV between 
atoms in nearest-neighbor sites. These facts hint that there is no bound state of H 2 in the bulk and/or that there 
are significant energy barriers that prevent the association of-dd to form H 2 in the bulk. The fact that Sieverts’ law is 
satisfied even at very high temperatures in certain materialsEI suggests that energy barriers are not solely responsible. 
On the other hand, high fluence hydrogen implantation experiments in graphite and some amorphous carbon filmsml 
suggest that H 2 formation may take place in the bulk after becoming saturated. This observation is not decisive since 
other interpretations cannot be ruled out, particularly in view of the high porosity of these materials. And, in any 
case, the pair interaction in the solid is undoubtably density dependent. At higher concentrations, the solid-hydrogen 


‘Tunneling of a hydrogen pair from bulk to vacuum is energetically possible in cases (i-ii) but not (iii). However, in case (ii) 
the decay time to a surface site is likely to be very short and the corresponding transition rate negligibly small. Differences in 
H 2 formation between cases (i) and (ii) may be important but are unlikely to be influenced by the bulk-vacuum tunneling rate. 
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systems tend, in some caseaEI, to undergo phase transitions to more ordered phases. We shall not consider such 
complications in the present analysis. 

Given the great degree of uncertainty we regard the binding energy of H 2 in the bulk and on the surface as free 
parameters. For illustration, let us consider the possibility that there exist bound ground states of H 2 (for a not too 
large separation) in or on the solid. We adopt the inequality {E c , E s } < discussed above for H. Furthermore, we 

consider only cases in which E b d > 0 and E[H 2 ] fc > 0 if H 2 levels exist in the bulk and similarly E s d > 0 and E[H 2 ] S > 0 
if H 2 levels exist on the surface. These inequalities imply that H 2 is energetically stable against dissociation in situ 
and that escape to the vacuum is energetically possible. 

If the condition that H 2 be bound in situ (with respect to H) is not satisfied, the existence of metastable bound 
states may still be possible in principle. The lifetime of such states depends on the characteristics of the corresponding 
energy barriers and is accounted for by the analysis presented in §|iy|. And, even if the condition is satisfied dissociation 
will still occur due to thermal activation and collisions with injected atoms, as discussed in greater detail in § |fv| . 

If the energy level of H 2 in the bulk and/or on the surface lies below that of H 2 in the vacuum, we may anticipate 
that the grain will become saturated with the molecules at sufficiently low temperatures. 

For each of the previously discussed cases (i-iv) one can identify all the possibilities consistent with the general 
assumptions. Only case (i) with 0 < E c < E s allows H 2 levels both in and on the grain. (Case [ii] allows H 2 levels in 
but not on the solid, case [iv] on but not in, and [iii] neither.) Figure |] illustrates one of three possible orderings for 
case (i). The thick arrows label different energetically allowed pathways for H 2 formation and release: H 2 formation 
in the bulk and on the surface, surface recombination, and evaporation of H 2 from the surface. Each reaction may 
involve an activation barrier with some characteristic height and width, as shown schematically in figure The 
association pathways (labeled 1 and 2) depend upon the H-H pair potential in the bulk and on the surface. The 
pair recombination from the surface (either quantum mechanical or thermal, labeled 3) is identical to the pathway 
discussed in the previous section. The evaporation of H 2 from the surface is labeled 4. 

We briefly review what is known about the association and dissociation kinematics of H 2 at the surface (i.e. paths 
2 and 3 and inverses). Direct computations of these processes are difficult and involve multi-.dipensional potential 
energy surfaces. They have been performed in a limited form for some hydrogen-metal systemstiZrES. Molecular beam 
experiments as well as self-consistent many-body calculations indicate .that for most metals there is an activation 
barrier for dissociation of width ~ 0.6 A about 2 A above-|the surfaceE£(E3. For simple and nobel metals the barrier 
ranges from about 0.2 eV (e.g., for Na) to more than 1 eVcJ. The repulsive interaction between the molecule and the 
surface is due to the molecule being closed-shell, and is similar to the repulsion found for closed-shell atoms such as He 
as they approach the surface. However, in the case of an H 2 molecule the antibonding state is shifted downwards and 
gradually filled as it approaches closer to the surface- thereby leading to a weakening of the attractive H-H interaction 
and the ultimate dissociation of the H 2 moleculecil. For transition metals (e.g, Ni, Pd) the activation barrier is 
appreciably smaller - about 0.05 to 0.1 eV - and for some systems (e.g., Ni(110)) the dissociation is non- activated. 
The reason for the small activation barrier observed in transition metals, as explained by Harris & Andersonc9, is 
that s electrons of the metal can occupy unfilled d states which are far more localized, --thereby reducing the Pauli 
repulsion between the H 2 core electrons and metal electrons. Self-consistent calculationso suggest that some metals 
may exhibit, in addition to the activation barrier of dissociation, a barrier for H 2 adsorption with a comparable 
height. Between the two barriers H 2 molecules can be trapped for a relatively long time, and this may account for 
the so-called molecular precursor state. For other materials there is no second barrier, but there is a small potential 
well just above the barrier-due to van der Waals forces that can give rise to physisorption of H 2 molecules at very 
low (< 20 K) temperaturesEJ. Large surface coverage seems to give rise to appreciably larger barriers, at least in the 
case of transition metals. For instance, the barrier for dissociation on a Ni(100) surface increases from 0.1 eV for a 
clean surface to about 0.6 eV for a surface with a full hydrogen monolayer coverageEI. Given these complications, we 
regard the activation barriers as parameters. 


E. Diffusion of Hydrogen in Solids 

Let the activation energies of diffusion of H and H 2 be Edh and Edh 2 . Most theoretical and experimental studies 
of H diffusion in solids have focused on diffusion in metals. We briefly review what is known about diffusion in solids. 
At sufficiently high temperatures, the dependence of the diffusion coefficient on temperature is well described by the 
Arrhenius law: 


D = D 0 exp(-E DH /kT), (4) 

where the preexponential factor, D a , is temperature independent _Jn the most rudimentary model of diffusion (for a 
review of the theory of diffusion of hydrogen in metals see e.g., ref.c3), the interstitial atom is supposed to be localized 
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at or about a given site. Diffusion occurs via a sequence of thermally activated jumps from one site to an adjacent 
one, in a random walk manner. The activation energy Ejjh is connected, in this model, to the height of the potential 
barrier. The jump frequency is expected to be of order of the zero point frequency for harmonic oscillations of the 
atom around its equilibrium position, denoted here by v a . Letting l denote the hopping length, the preexponential 
factor can be approximated as 


D a ~ 10 3 <7(£/A) 2 (V o /10 13 s 3 ) cm 2 s % (5) 

where q is a geometrical factor which depends on the lattice type. 

The values ofpQ 0 and Edh have been measured for some materials using various experimental techniques (for 
a review see refJl3). For most metals, as well as for some nonmetals for which there are experimental data, D a 
lies in the range 10~ 4 to a few times 10~ 2 cm 2 s _1 , in accordance with Eq. The activation energy appears 

to vary appreciably from one substance to another, ranging from about 0.043 eV for vanadiumES to ~ 0.5 eV for 
titaniumEj. It should be noted that most of those measurements have been taken at relatively high temperatures and 
in a fairly narrow range. Moreover, for some materials (e.g., iron) there are quite large uncertainties, particularly at 
low temperatures. The determination of hydrogen diffusivity in carbon is difficult because of the presence of deep 
trapping sites that appear to control the mobility at low hydrogen concentrations; the diffusion energy inferred from 
experiments using low concentrations of hydrogen is about 4 eV, which is comparable with the trap energy. The best 
estimate of hydrogen diffusivity in graphite (when trapping sites are excluded) according to Causeycil is Edh — 2.8 
eV, D a ~ 0.93 cm 2 s _1 . 

At low temperatures (~ 200 K for metals) quantum effects play a role, and deviations from the Arrhenius law 
are expected, and in some cases have been observed. (1) At extremely low temperatures, the interstitial eigenstates 
should fulfill Bloch’s theorem and form a band. The diffusion process is then band propagation similar to electron 
conduction. The rate, limited by interaction with phonons and lattice! defects, decreases as temperature rises. No 
indications of band propagation for hydrogen have been observed yetEHI. (2) At somewhat higher temperatures, the 
decay rate of the band states increases, and ultimately approaches the band width. As a consequence, at temperatures 
larger than some critical temperature the interstitial will be localized about a specific site, as discussed above. In this 
temperature regime, overbarrier transitions in the manner described above are negligible, and thermally activated 
tunneling from one site to another may dominate the diffusion process. In the small polaron theoryc3 one treats 
the tunneling matrix element as a small perturbation in the total Hamiltonian, it is then possible to compute the 
transition rate (i.e., hopping frequency) using time-dependent perturbation theorylkx The results of such calculations 
show that well below the Debye temperature the transition rate should satisfy a T 7 law, which has not yet been verified 
experimentally. At temperatures greater than the Debye temperature, the transition rate obeys the Arrehnius law 
with an activation energy that corresponds to the energy difference between an occupied site and a vacant site (note 
that the slope of the log D vs. 1/T curve is different than that in the high temperature regime where overbarrier 
transitions dominate) and a preexponential factor that involves the tunneling matrix element (and therefore should 
be strongly isotope dependent). 

From the foregoing discussion it is expected that at sufficiently low temperatures the diffusion coefficient will 
generally have non-exponential temperature dependence. The temperature at which this happens depends on the 
width of the potential barrier and other details of the interaction of the H atom with the lattice. Unfortunately, 
measurements of the diffusion coefficient at low temperatures are extremely difficult, and the value of D below about 
200 K is poorly known. 

The mobility of hydrogen on surfaces is some 10 - 15 orders of magnitudes larger than that in the bulk by virtue 
of the small barriers between surface sites. Recent measurements indicate that there is indeed a sharp transition for 
surface diffusion from a high temperature regime, wherein D obeys the Arrehnius law, to a low temperature regime, 
where Lbis essentially temperature independent. For example for W(110) at zero coverage this transition occurs at 
~ 150 Kcj. However, there is evidence for a strong anomalous isotope effect, as well as some other complications, 
which are not well understood at present. 

At high concentrations, hydrogen diffusion in the solid may be altered significantly, owing to the increasing strength 
of the self interaction between the hydrogen atoms. 

In view of the large uncertainty in D referred to above, we shall treat the mean hopping time between bulk sites, 
defined as th = ^ 2 /D, as a free parameter. 


IV. A MODEL OF MOLECULAR HYDROGEN FORMATION IN GRAINS 
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A. Basic Equations 


We distinguish sites according to whether they lie in the bulk or on the surface of the solid and whether they are 
elements of the regular solid (“clean” sites) or correspond to centers of altered binding (impurities, vacancies, and 
defects collectively called “impurities”). In our model a site may be empty (</>), occupied by a single H with energy 
£J[H] or occupied by H 2 with energy we do not allow multiple occupancy or excited energy levels. Denote 

by Paa(x,t) the number density of bulk sites of type a ( a refers to clean bulk sites or to impurity sites of type k 
denoted Ik) that are occupied by species a (a = H, H 2 ) at position x at time t. Similarly, denote by a aa (x,t) the 
corresponding areal number density of surface sites. 

We assume that the average number of H and/or H 2 per grain is much larger than 1 so that rates of relevant binary 
reactions can be meaningfully expressed in terms of the average number densities of the reactants. 

We assume that species diffuse in the bulk and along the surface of the grain and that exchange occurs between 
the gas phase, the solid bulk and the surface. We express the rate of change of the number densities by the following 
coupled set of equations: 


^ pact 

~df 


^ {H act^ P act) Pact ■> 


dcr aa 

dt 


^s{I^aa^ s ^ aa ) d aQ ., 


( 6 ) 


where D aa and D s aa are the bulk and surface diffusion coefficients, respectively, and V s denotes gradient along the 
surface. The terms on the right hand side give the rate of change of density from external sources (we denote for 
future use the injection from the gas phase of species a into sites of type a in the bulk and on the surface by s aa 
and s aa , respectively) and from internal rearrangements of site occupancies. Internal changes in the bulk number 
densities occur from (i) chemical reactions taking place between bulk species or between bulk and surface species, 
(ii) particle exchanges between bulk and surface sites, and (iii) particle losses to the vacuum (either via evaporation 
or recombination). Rate coefficients for these processes depend on the parameters of the grain (temperature, energy 
levels, and so forth, including occupation probability). 

Let S a = (1,2) for a = (H, H 2 ). Summing eqs. (||) over the entire grain, over all species and filled sites yields the 
total rate of change of hydrogen nuclei 


d 

dt 





dV + 



-‘It- 1 

L in u H—out H-2—out 


(7) 


where f” 1 is the rate for atoms to strike the grain (eq. |]]), t~j^_ out is the loss rate of atoms and t^_ out is the loss rate 
of molecules. The loss terms include loss by evaporation (thermally and quantum mechanically), loss induced by fast 
collisions (direct ejections, pickup reactions, enhanced recombinations) and loss by pair recombination (thermally and 
quantum mechanically). In steady-state the flux of atoms intercepted by the grain is balanced by the total loss rate of 
atoms and molecules. Efficient molecule formation requires the total loss rate of molecules by all channels to exceed 
that of atoms. The relevant rates depend on the occupation numbers of the different species, which are determined by 
the various physical processes taking place in and on the grain. Below we explore the solutions to the above equations 
in different regimes of the parameter space, and elucidate the conditions under which effective molecular formation 
may take place. 

In the examples presented below, we further assume that the occupation numbers are homogeneous within the 
bulk and homogeneous on the surface and that the grain temperature T g is uniform throughout. In other words, we 
describe the chemical distribution in terms of two zones (bulk and surface). This approximation is justified whenever 
the diffusion timescale is short compared to all other characteristic timescales (for example, the injection time) and 
may also be reasonable under less restrictive conditions (for example, if the injection process is homogeneous). We 
define the total number of grain sites on the surface and in the bulk (occupied a or empty <f>) 


N= E 

ol,ol,4> 


PactdV + 



( 8 ) 


the fraction of A f that are in the bulk, of type a, and occupied by species a 




the fraction of A f that are on the surface, of type a, and occupied by species a 
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aa 


= U ~ 1 


1 




and the bulk and surface injection rates 


S 


aa 



S a 



(10) 


( 11 ) 

( 12 ) 


Steady state solutions are obtained by numerically integrating the time dependent rate equations (jlj) (integrated 
over the appropriate zone) until equilibrium is achieved. For the specific contributions we have included forward and 
backward reactions related by detailed balance for each process that involves thermalized surface and bulk reactants. 
For reactions driven by impinging fast particles and for particle loss to the vacuum we have included only the pathway 
of interest. We have not included any direct bulk to vacuum loss rates even when quantum tunneling might be 
possible. 

Our analysis becomes inapplicable when the mean residence time of an hydrogen atom in the grain, At becomes much 
shorter than ti„.. In this case the probability of finding two atoms simultaneously inside the grain is approximately 
A t/ti n , assuming that injection of atoms is a Poisson process. Consequently, even if the rate of molecular formation 
is shorter than At, at most a fraction A t/ti n of injected atoms will be converted into molecules. 


B. H2 Formation in a Clean Grain 

We consider the simplest case, namely molecular formation in a grain free of impurities ( ni a = 9j a = 0). We begin 
by giving explicit expressions for the rate coefficients in terms of the various parameters involved. 


1. Reaction Rates 

Below we list the processes included in the model. As a general strategy, we have parameterized the rates in terms 
of explicit functions of nn, njy 2 , 9h and 9h 2 (“fractional occupancies”) times characteristic rate coefficients (of the 
forms t , aS). These latter depend primarily on T g and various grain energy scales but also on occupancies when 
saturation is approached. (Later results are explicit and most useful when this implicit dependence on fractional 
occupancy is absent or weak but remain correct implicit solutions in all circumstances.) 

Let g s and gb be the fractions of all sites that are surface and bulk sites respectively (i.e. ~ A/ --1 / 3 and 1 — A/" -1 / 3 , 
up to a geometrical factor; g s + gb = 1). The fractional occupancies satisfy 0 < nn + tih 2 < gb and 0 < 9h + 9h 2 < gs- 
We denote an empty bulk site b(f>; the probability that a bulk site is empty is B b = 1 — (nn + nti 2 )/ 9b- Likewise, we 
denote an empty surface site scj>; the probability that a surface site is empty is B s = 1 — (0 H + 9n 2 )/gs- 

The reactions are listed and the forms for the characteristic rate coefficient are given immediately to the right, 
followed by a short description of each process and the expression for the rate of change (always expressed as per site 
for the J\f grain sites). 

• Particle Exchange 


ba + s<f> —*■ b<j) + sa t bs l a 

sa + 6</> —► scj) + ba t~ ba 

— Thermally activated diffusion of a species a = H or H 2 from the bulk to the surface and vice versa. The 
fluxes are proportional to the product of bulk (surface) concentration and the number of vacant neighbor 
surface (bulk) sites. The rates per site may be expressed as t b ^ a n a (g s B s /gb), and t~ ba 9 a B b where tb sa and 
t s b a are characteristic timescales to hop to nearest neighbor sites. The hopping timescales depend on the 
grain temperature, the characteristic binding energies and the form of the interaction. We assume exchanges 
between neighboring lattice sites; order unity range variations are accommodated in the definition of the 
characteristic timescale. Detailed balance implies t~ ba = exp {(E s [a\—E b [a\)/kT g }t b ^ a . If E b [a\ — E s [a] =0 
then tbsa — t s ba ^ tha- 
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H* + 6H + s<j> —► H*' + b(j) + sH 
H* + sH + b<j) —> H*' + scf) + 6H 


0!bs B 

a sb S 


— Collisional displacement of H atoms from the bulk to the surface by energetic incident atoms (H*) and vice 
versa. The rate coefficients are taken to be proportional to the total injection rate per site, S = 1/Aft in , 
with efficiencies which are denoted by a bs and a sb , respectively. The rate per site for transport from bulk 
to surface takes the form a bs Sn,H(g s B a /g b ) and from surface to bulk, a sb S9HB b . 

The efficiencies depend only on the energy of the incident particle and the energy level of the various sites 
as long as the density is smaller than the inverse volume sampled by the incident particle. Once the grain 
becomes saturated the efficiency of collisional displacement will depend on density, for example, when one 
dislodged atom can dislodge additional atoms. 

A plausible limiting efficiency may be derived as follows. If in slowing down, an incident atom visited every 
site in the grain once, then each occupying atom suffers ~ 1 displacement. A surface atom is moved into 
the bulk with order unity probability; a bulk atom is moved into the surface with probability g s . From 
the form of the above rate expressions, this corresponds to a bs and a sb ~ TV". Thus, we assume that the 
efficiencies range from zero to roughly TV". 

• H 2 formation and destruction 


6H + 6H - 

-> &H 2 + b<j) 

t~ l 

L Fb 

sH + sH - 

—> sH 2 + s(j) 

t~Fs 

6H + sH - 

-> &H 2 + scj) 

1 

L Fsb 

6H + sH - 

-> sH 2 + b(j) 

r 1 

L Fbs 


— Molecular hydrogen formation in the bulk and on the surface (reaction pathways labeled 1 and 2 in figure j|). 
We define t-Fb (tp s ) to be the H 2 formation time, given that two H atoms are in neighboring, bulk (surface) 
sites. We define tFbs ( tFsb ) to be the time of H 2 formation on the surface (bulk) by recombination of 
subsurface and surface atoms. If H 2 formation in or on the grain is forbidden then tFi = tFij = 00, while in 
the absence of an activation barrier tF b ~ th- Generally tFi are finite and depend on the height and width 
of the corresponding activation barriers for H2 formation in the solid, as depicted in figure [zj. The rate 
per site for H 2 formation in the bulk is [tplnn + tplb^HjiT-H / g b , and that for H 2 formation on the surface 
\Pf s {Qh / Qs) +tFbs( n H/9b)\0 H - As above, we have assumed that the hydrogen pair interaction range is one 
lattice spacing but the variation may be absorbed into the definition of tpi- 

H* + 6H —» 6H 2 a Fb S 

— H 2 formation in the bulk by direct recombination of injected atoms with interstitial atoms. The rate of 
change of nn 2 is given by ctFbSriH where the efficiency aFb depends only on the energy of the incident 
atom and the site binding energy. Again, aF b ranges between zero and TV". If H2 formation in the bulk is 
forbidden, then aFb = 0. 


H* + sH 


sH 2 

uH 2 


C^FsBSfs 
aF S S(l - Sfs ) 


— H 2 formation on the surface by direct recombination with the incident atom. The corresponding rate of 
change of 9p 2 is given by aFsB9p where the efficiency is Qjf s . A fraction Sf s of the formed molecules 
are assumed to remain on the surface, while the rest are immediately ejected from the grain. Following 
DuleyHS, we shall refer to the latter process as prompt reaction. 


6H 2 + b(j) — 

611 + 611 

l Db 

sH 2 + S(j) — 

-> sH + sH 

r' 

Dbs 

6H 2 + s(j) — 

-► 6H + sH 

sH 2 +bcj) — 

-► 6H + sH 

L Dsb 


— H 2 dissociation in or on the solid owing to thermal activation and quantum tunneling (in the case where 
metastable states exist). Detailed balance implies that the ratio of formation to dissociation rate tjq ftp] = 
exp(2 E l d /kT g ), i = b,s. We take the dissociation rate per site to be tJ ) \nH 2 B b in the bulk, t Ds e H 2 B s on 
the surface. Likewise the rate coefficient for dissociation of H 2 in the bulk giving surface and subsurface H 
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is ^Dbs = ^Fsb ex p{ — (2-E 1 ^ + E[H] S — E[H] b )/kT g }; for H 2 on the surface the corresponding rate coefficient 
is = ^~Fbs exp {-{2E s d + E[H] b - E[H] s )/kT g }. The rates become t~l s nH 2 (g s B s /g b ) and tpl b e H2 B b , 
respectively. 


H* + &H 2 + bcj) - 

-> H* +6H + 6H 

PbS 

H* + 5 H 2 + scj) — 

-> H*' + sH + sH 

&(1 ~q)S 

H ¥ + SH 2 + scj) — 

-> H*' +vE + vR 

P s qS 


— H 2 dissociation in or on the solid by injected atoms. The bulk and surface dissociation rates per site are 
PbSnH 2 B b and /3 s S9h 2 B s where (3i gives the average number of dissociations per target molecule per site 
as the atom slows down for i = b,s. Further, we suppose that a fraction q of H 2 dissociations on the surface 
results in the immediate loss of the H atoms thereby produced from the grain. 

• Particle Losses 

— Thermal evaporation of H atoms from the grain’s surface at a rate 7 h^h- 
— Collisional ejection of H atoms from the grain’s surface at a rate q: sv hS9h- 

— Evaporation of H 2 molecules from the surface (thermal and quantum traversal of pathway 4 in figure [^) at 
a rate given by r yp 2 9h 2 ■ 

— Collisional ejection of H 2 molecules from the grain’s surface at a rate a sv n 2 S9H 2 - 

— Pair recombination of adjacent surface atoms (thermal and quantum association via pathway 3 in figure ||). 
with a rate given by g~ 1 T 2 h^h- We henceforth refer to this process as “pair evaporation.” 

For our models we need only the sum of various individual processes. These are: 

• The rate coefficient for bulk to surface exchange is = t bsH + ctb s S ; the total rate is nH{g s B s / gb)r b ^ H . 
Likewise, the rate coefficient for surface to bulk exchange is r~ bH = t~ bH + a s bS; the total rate is 9nB b T~ bH . 
The rate coefficients for H 2 exchange are completely analogous. 

• The dissociation rate coefficient on the surface is Tp s = t~f^ s + (3 s S] in the bulk r ^ 6 = tp b + Pb‘S', and for surface 
plus subsurface products t~^ bs and t~^ sb . We also define the part of the surface dissociation rate that leads to an 
increase of 9h (as opposed to immediate escape) TDs T Ds QPsS. 

• The rate coefficient for thermal and collisional loss of H from the surface (all linear in 9h ) is F p = 7 h + a sv nS. 
The rate coefficient for H 2 loss by thermal, quantum and collisional loss is Th 2 = 7 h 2 + oi sv h 2 S- 

• The H 2 formation rate in the bulk is riH{tplnH/gb + tpl b 9n/gb + cuFbS ); the thermally driven H 2 formation 
rate on the surface is 9H{t~p l s 9n / g s J ri~p\ s n H/gb ); the H 2 formation rate by collisions and by pair evaporation is 
9h{chfsS + gj 1 r 2 _f/0_f/). 


2. Relation of Rate Coefficients to Energy Scales 

We now turn to the evaporation rates. We suppose that the interaction potential of a single hydrogen atom with 
the surface is similar to that shown schematically in figure |lj. with E e > 0. The evaporation rate of an H atom from 
a surface site is 

7 h = v 0 exp{—i?e/fcTg}, (13) 

where v Q ~ 10 13 s _1 is the characteristic frequency of oscillation of an H atom about its equilibrium position. 

The pair evaporation rate, T 2 h, depends on the shape of the pair surface potential. Atom pairS|-&an escape from 
the surface by either hopping above the finite activation barrier or if E c > 0 by tunneling through itll3. When E c > 0 
recombination can take place even at zero temperature and tunneling dominates the overbarrier rate at sufficiently 
low temperatures. In the following, we assume that the initial separation between the recombining atoms is of the 
order of the H-H interaction range (about 1-2 surface sites). 

We consider first a grain with E c > 0 at zero temperature, as shown in figure ^](i) and (iv). Let R be the vector 
position of the center of mass of the pair system, r one-half the relative position vector of the nuclei, and V (r, R) the 
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corresponding adiabatic potential. To zeroth order, the penetration factor of an atom pair of energy E along a given 
path joining the points a and b in this six-dimensional space, is given in the WKB approximation by 


G s = 2 ( 2 m 2 H /n 2 ) 1 / 2 f dS[V (r(5), R(5)) — E } 1 ^ 2 , 

J a 


(14) 


where dS is a line element, and V(a) = V(b) = E. The transmission coefficient can be determined by summing e~ Gs 
over all paths joining the points a and b , and the corresponding tunneling rate is obtained by multiplying the result by 
v 0 , the attempt frequency. Lagrange’s equations describe the path of minimum penetration factor for fixed endpoints 
a and b. The dominant contribution to the tunneling rate comes from that path and to a good approximation at zero 
temperature the surface recombination rate is 

V 2H ~ v 0 e~ G , (15) 


where G denotes the appropriate value of Gg. To simplify the analysis further, we consider only paths along which r is 
constrained to be parallel to the surface. We assume that the potential energy surfaces are co-planer for fixed separation 
so that the potential depends only on the two coordinates z and r, where 2 is the height above the surface (i.e., the 
component of R normal to the surface) and r = |r| is one-half the distance between the hydrogen nuclei. Then the 
penetration problem reduces to tunneling through an effective one-dimensional potential v(z) = [V (z)—E] [1 +{dr/dz) 2 ] 
where ( dr/dz) describes the appropriate path. In view of the uncertainties in the 2D potential energy surfaces, we 
approximate v{z) by an inverted harmonic oscillator potential with barrier height AE and width z a . Then 

G h = 48.5(A E/eV) 1 / 2 {z 0 /A) (16) 

and eq. (^5|) gives the zero temperature recombination rate. The result is easily generalized to finite temperatures. 
Averaging the tunneling rate over the Boltzmann distribution while regarding the transmission coefficient as unity 
for E > AE yields 


F 2 h 


Vo 

G h - A E/kTg 


G H e~ AE,kTg 


(A E/kT g )e~ GH 


(17) 


The rate is strongly enhanced when the temperature approaches the tunneling temperature, T t = AE/kGn — 
2.4 x 10 2 (A E/eV) 1 / 2 (z 0 /A) -1 K, and approaches the classical limit above Tt- 

For materials with E c > 0, the pair evaporation rate is given by eq. JTt| ) with AE = E r (see figure [|[i,iv]), and Gh 
the penetration factor at energy E c , i.e. the ground state energy of the trapped particles. The limiting behavior is 


2 H 


-Gh min(2i-,l) „ n 

= v a e V J g ) it E c > 0. 


(18) 


For materials E c < 0, thermal activation is required (see figure j^[ii,iii]). Only atom pairs with total energy 
> E[H 2 ] v can tunnel through the barrier and we find that the surface recombination rate is eq. (0 reduced by a 
factor exp (—E c /kT g ). Here, AE = E a (see figure |l][ii,iii]) and the interpretation of Gh is that it is the penetration 
factor of an atom pair having energy E[H 2 ] V . The limiting behavior is 


F 2 h = v 0 


Eg 

O kTn (_ 


-G„min(%,l) i£E _ 


(19) 


In the case in which H 2 formation in or on the solid takes place, the evaporation of H 2 molecules from the grain 
surface proceeds in an analogous fashion. If E[H 2 ] S > E[H 2 ] v , as assumed for figure then H 2 escapes by both 
overbarrier and underbarrier transitions. On the other hand, if E[H 2 ] S < E[H 2 ) V then H 2 requires thermal activation 
to escape. Like T 2 h, the H 2 evaporation rate, 7 # 2 , depends on the energy barrier and whether or not tunneling can 
occur at zero temperature. 


3. Rate Equations 

We suppose that the occupation numbers are homogeneous in the bulk and (separately) on the surface. As already 
mentioned above, this assumption is justified when the diffusion time is short enough. The equations governing nn, 
n H2 , 0h 2 and Oh are given respectively by, 

n H = S H - [t^ s h{9sB s / g b ) + 2 a Fb S]n H - ^9 b l tp\n 2 H + T~ b i H B b 0H 

+ t Dbs9sB S /g b )nH 2 +t Dsb B 0 h 2 - g b (^Fsb d” ^Fbs) nH ^H, (20) 
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( 21 ) 


nH 2 = -[ T Db Bb + {a sB s / g b )t D \ s + T bs \j 2 (g s B s / g b )}n H2 + a Fb Sn H 

+ t Fb9 b ln H + T sbH 2 B b ®H 2 + g b 1 t F 1 sb riH9H, 


9h 2 = T bs H 2 nH 2 {g s B s /g b ) +t F \g s 1 9 2 H + a Fs S Fs S9 H + g b 1 t F l s n H 9 H 

~ ( T sbH 2 B + + t D s B +t Dsb B )9h 2 , 

and 


( 22 ) 


9h = S h + r bs 1 H n H {g s B s /g b ) - (r ; 1 
L Fbs )g b 1 n H 9 H + (2 t Ds 


~ (tFsb + 1 


-1 


sbH 

-1 


B b + 2 a Fs S + Y h )9h - 2(r 2 h + t F ])g, 
B +t Dsb B )9 h 2 + (g s B s /g b )t Dbs n H2 - 


rjl 

y H 


In steady state fin = hH 2 = 9h 2 = 0 F = 0. Then the sum of eqs. (EG) and (EG) and twice eqs. 


and 


(23) 

yields, 


S = S H + S H = t H _ out + 2 t H2 _ out 


where the loss rate of H to the vacuum is 


L H—out 


= Yh9h + 2 B s (t D l — t d I)9h 2 = Y h 9h + 2 q(3 s SB s 9 H21 


and the loss rate of H 2 to the vacuum is 

^ 2 -out = Fh 2 0h 2 + g~ 1 T 2 H9 2 H + a Fs S0 H (1 - S Fs ). 

The efficiency of H 2 formation is the fraction of hydrogen released from the grain in molecular form: 


£ = 


2 1 


-1 

H2~OUt 


t 


H—out 


2 1 


-1 

H^—out 


= 1 - 


b H—out 

S ' 


(24) 


Let us examine the solution to the above equations in various limits. First, suppose that the timescales for 
H 2 formation and thermal dissociation on the surface, t Fs and td s , are much shorter than any other timescale 
involved. We then anticipate the H and H 2 to be in equilibrium on the surface. Indeed, in the limit where all 
the rates except the formation and dissociation rates tend to zero, equations (|22|) and ( pj| ) yield, B s (9h 2 /g s ) = 
( T Ds/^Fs)(^/ffs) 2 = e 2 ' B ^/ fcT s(0#/g s ) 2 , which is just the appropriate mass action law (note E s d is the dissociation 
energy per atom). Similarly, when H 2 formation and dissociation are the dominant processes in the bulk, eqs. ( |2C| ) 
and (21) yield, {uh 2 / gb)B b = (TDb/tFb)(nH / gb) 2 = e 2Ed ^ kTg (uh /gb) 2 ■ Next, suppose that the surface (bulk) to 


bulk (surface) diffusion times are very short compared with any other timescale. Upon taking the limit where all 
rates except the exchange rates tend to zero in eqs. (^) and (||^), we recover the equilibrium distribution of H: 
B b (,9 H /g s ) = B s {t sbH /t bsH ){riH/gb) = B s exp{(E[H] b - E[H] s )/kT g }(n H /gb)- Finally, in the limit where the H 
evaporation time is much shorter than quantum pair recombination, H 2 formation and dissociation on the surface, 
and transfer times from bulk to surface and vice versa, we obtain from eq. (|23|), Oh = Sh/Yh , as expected. Likewise, 
when the pair recombination rate dominates, eq. (23) gives 9h = (gs$H/^Y 2 h) 1 ^ 2 ■ 

In the following, we shall neglect the densities in the statistical factors, as well as the surface-subsurface reactions. 
Specifically, we set B s = B b = 1, t F ] b = t F l s = t~^ sb = t~j^ bs = 0. This approximation enables us to solve eqs. ( |2(i| ) 
through (|2^) analytically in various regimes. The neglect of the densities in the statistical factors is certainly justified 
when the grain is unsaturated (i.e., UH+nH 2 « g b , 9h +9h 2 « gs), & condition we check in our solutions. The full 
numerical solution of eqs. (|2l])-(|23|) indicates that this is a good approximation even for relatively large concentration 
numbers. In ignoring cross terms between the surface and bulk we are essentially assuming that the mixed rates are 
not larger than the corresponding rates in both the surface and the bulk. We have checked our results in a variety of 
limits and find no qualitative difference for the parameter ranges we have explored. 


4- Equations without H 2 Formation/Destruction in the Bulk 

We now explore the conditions required for effective H 2 formation in different regimes. First, consider a grain for 
which neither formation nor destruction of H 2 is possible in the bulk. This corresponds to the limit a Fb = t F l = 
t~D b = (3 b = 0 in eqs. (|2(i| ) - (f23[). Then eq. ( ^0|) can be written in the form, 
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eq- (H, 


n H (ai, a 2 , 03 ,a 4 ) = a 3 0 ff (ai,a 2 ) + < 24 , 


(25) 


n-H 2 ( a l, a 2, « 5 , a 6 , “7) — a 7 ^ff 2 ( a l, ° 2 , a 5 ; a 6 ), 


(26) 


and eq. (H), 


The sum of eqs. 


Here 


6> ff2 (ai,a2,a 5 ,a 6 ) = a 5 0^(ai, a 2 ) + a 6 0j?(ai, a 2 ). (27) 

and (H|) gives the following dimensionless equation for the surface occupation probability: 

1 - a 1 e H (a 1 ,a 2 ) - a 2 6 2 H (a 1 ,a 2 ) = 0. (28) 


ai = ThS 1 + 20 :^ 


1 , 8fsT Ds 

OtFs = a.Fs 1 — 


T n.s + Tff, 


«2 — 2 g s 5 (r 2 //+t i?s ); 

°3 = {%/ 9 s){T bsH /T sbH ), 
°4 = (W g s )T b sHSH, 

9s ^Fs 

05 = —ZI- 

T; 


f -1 — +- 1 I 1 _ 

''P’s ~ Ts \ 1 -1 




T Ds‘ + r ^2 


Ds T 
_ aFs^FsS 

° 6 _ _-i 1 r ’ 
t Ds + 1 #2 

°7 = (gb/gs)(r b sH 2 /T sb H 2 ), 


(29) 


are dimensionless parameters that depend on the basic physical parameters as stated explicitly above. 

As seen from eqs. ( |28| ) and (^fj), the surface occupation number, Oh, depends on only two of the seven dimensionless 
parameters of the problem. For example, only the total injection rate (S') but not the individual bulk ( Sh ) and surface 
( Sh ) injection rates are relevant for the value of Oh- Since there is no interconversion of H and H 2 inside the grain, 
in steady state every atom intercepted by the grain must eventually escape through the surface. We will exploit the 
dependence of Oh on ai and a 2 in the discussion in the next two sections. Note that in the simplest case with no H 2 
formation on the surface (tp s = 0 ), no prompt reaction (aj? s = 0 ) and pair evaporation only by quantum tunneling we 
have — lnai oc — lnT# oc E e /kT g and — lna 2 oc — lnr 2 # oc Gh- The essential microphysical parameters that govern 
the surface H density are E e and Gh- We will begin by characterizing the solutions to the steady-state problem 
directly in terms of possible combinations of E e and Gh ■ We will then argue that more general cases are directly and 
simply described in terms of similar regions expressed as functions of — lnai and — lna 2 , the analogues of E e and 
Gh- 

Although the surface density depends solely on a± and a 2 , all the dimensionless parameters are needed to find a full 
solution. That is important because we must check for saturation of surface and bulk components as well as calculate 

e. 


5. H 2 Formation by Pair Evaporation and/or Prompt Reaction 

To elucidate the role of pair evaporation and prompt reaction, we further assume a grain for which H 2 is unstable 
on the surface. This corresponds to the limit tf> s 0, ths/tds —> 1 in eq. (HI)- Then, a 2 = 2g/ 1 T 2 H/S and 
a -5 = a -6 = 0. Eqs. ( |2(i| ) and (H|) then implies, nn 2 = 0h 2 = 0. Now, when 0,4 > g b the injection rate of H atoms 
into bulk sites, Sh, exceeds the rate at which hydrogen is released from the bulk, so the bulk will become 

saturated regardless of the values of the other parameters, as may be seen from eq. (|25|). If a 4 << g b bulk occupancy 
is controlled essentially by H transfer from the surface to the bulk rather than injection of H in the bulk, that is, by 
the first term rather than the second on the right hand side of eq. (|2E|). The grain will remain unsaturated if the 
evaporation, surface recombination or prompt reaction rates are sufficiently large. (The efficiency £ = 1 — ThOh/S 
(see eq. |24|]) is also independent of 04 .) Provided we are not interested in the exact value of nn so long as it remains 
below g b we can proceed by setting <24 = 0 in eq. (HI)- The solution for nn is then exact in situations in which 
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the ambient gas is cold, such that H atoms can only stick to the surface (i.e., Sh = 0, Sh = S), and is a good 
approximation in general, provided that the rate at which H atoms are intercepted by the bulk is much smaller than 
the characteristic rate of transfer from the bulk to the surface. 

When 0,3 < gb/g s , or equivalently TbsH/TsbH < 1, the steady-state distribution of hydrogen in the grain always 
satisfies nn/gb < ^n/gs- Consequently, when the surface reaches saturation the bulk remains unsaturated. When 
03 > gb/gs-, bulk saturation always precedes surface saturation. Let us define a critical surface concentration above 
which the grain (either bulk or surface) is saturated: 


J 1 ; TbsH/T s bH< 1 surface saturated 

\ T s bH/nsH ; TbsH/TsbH > 1 bulk saturated. 


(30) 


If thermal diffusion alone is responsible for surface-bulk fluxes, then TbsH/TsbH = hsH/t s bH = e~ < ' Es ^ Ea ^ kTg and 
the top line applies for E s > E c , the bottom for E s < E c . The critical surface concentration has the simplest form 
(9 sat = 3 s) when E s > E c . The condition for either bulk or surface to be saturated is temperature dependent when 
E s < E c or whenever the rate coefficient includes collisionally driven exchange. At several points in this paper we 
will make the further assumption that the equations governing surface concentrations continue to apply even after 
the bulk has become saturated. Formally, the blocking factors that we have dropped would ensure mathematically 
well-defined solutions. The critical assumption is really that the basic grain parameters, measured or calculated for 
low occupancy continue to remain accurate in the saturated limit. 

On substituting 9 sa t into eq. (^), we obtain the equation of a surface in the space of physical parameters that 
separates the regimes of saturated and unsaturated grain. Specifically, the condition that the grain will be unsaturated 
reads, 


^2 h > f _9sS_\ 

V a \ 29 S atVo J 



2 aFs 



(31) 


with 


— = exp (-E e /kT g ) + (32) 

Vo V 0 

where eq. (^) together with the definitions of the various physical parameters given above have been used. This 
inequality may be viewed as a relation between r 2 # and E e /kT g when other parameters are fixed. In the absence 
of collisional ejections and prompt reactions the grain will be unsaturated either if T 2 H > gsS/29^ at or if E e /kT g < 
In (y o 0 sa t/S). 

In a similar manner, we can derive an equation describing a family of isoefficiency surfaces, where the efficiency is 
defined in eq. (pi|): 


Y2H 

v 0 


g s Yn 

2(1 -£) 2 Sv a 


(£Th 


2(1 -£)a Fs S). 


(33) 


Note that the isosurfaces are independent of the timescales T s bH and Tb S H (and 9 sa t), as one might expect, since (i) 
evaporation and H 2 formation processes are restricted to the surface in the example presented here, and (ii) as already 
emphasized, the surface concentration depends on the total injection rate alone. In the absence of collisional ejections 
and prompt reactions eq. (^3|) assumes the simple form, — \n{Y 2 H/v 0 ) = 2E e /kT g + ln[2S(l — £) 2 / (£g s v 0 )\. 

Projections of the surfaces governed by eqs. ([H]) (when the inequality sign is replaced by the equal sign) and ( fi^ j 
onto the plane defined by — \n(Y 2 H/v 0 ) and E e /kT g are exhibited graphically in figure |], for g s = 10 -2 , and different 
choices of the remaining parameters (note that the dependence of the curves on g s as well as the other physical 
parameters besides — ln(r 2 fl-/z/ 0 ) and E e /kT g is logarithmic). The regimes of saturated grain, effective and ineffective 
H 2 formation are indicated. The regions indicated as effective and ineffective in the figures are separated by the 
curves along which the efficiency £ = 1/2. In the regimes indicated as effective the pair evaporation rate exceeds 
the evaporation rate of H and in those indicated as ineffective the opposite is true. Windows (a) and (b) present a 
system not subject to direct collisional processes ( a SV H = = 0). The region enclosed by the solid line in case (a) 

corresponds to a saturated surface, whereas the regions enclosed by the curves labeled by values of Tbsu/^shH > 1 
correspond to bulk saturation. Note that for any given substance, 9 sat is variable, as discussed above. (In a later 
section we speculate upon the steady state of grains saturated with hydrogen, and discuss several plausible scenarios 
for H 2 formation in such systems.) 

It should be noted that T 2 # is constant only in the limit of zero temperature quantum tunneling. The pair 
evaporation rate generally depends on the grain temperature and T 2 # —> v 0 as T g —► 00 . For E c > 0, the limiting form 
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for pair recombination (eq. |Q) implies — ln(r 2 7 f/V 0 ) — (E r / E e )(E e /kT g ) at grain temperatures above the tunneling 
temperature T t , and — In(r 2 j 7 /V 0 ) — Gh below; for E c < 0 the limiting form (eq. [^9|) implies — ln(r 2 .H-/i/ 0 ) — 
[( E a — E c )/E e \(E e /kT g ) above, and — ln(r 2 if/i/ 0 ) — Gh — E c /kT g below. In figure ||, we show the dependence of pair 
evaporation rate on grain temperature (—lnT 2 h/vo is plotted against E e /kT g ) for two materials. The dashed line 
describes a typical transition metal with E e = 2.55 eV, E c = —0.3 eV, E a = 0.2 eV, and Gh > 7. The dotted-dashed 
line describes carbon (with E c > 0) using the same parameters adopted in §[v| namely E e = 0.3 eV, E r = 0.4 eV, 
and Gh = 30. The solid lines for efficiency and surface saturation were computed using the same parameters as in 
figure ||(a) with t bsH/ t s i,h < 1. Note that the conditions for bulk saturation are not displayed in this figure. In the 
case of carbon we find that bulk saturation occurs (for E s = —0.2 eV, S ~ 10~ 13 s _1 , and the above choice of the 
remaining parameters, see §0 for a detailed discussion) when E e /kT g > 7. For both of these putative materials the 
surfaces remain unsaturated. The threshold temperature below which effective H 2 formation takes place is ~ 10 3 K 
for the transition metal and ~ 90 K for carbon. 

The effect of collisional ejection of H from the surface and prompt reaction is illustrated in figure windows (c) 
and (d). The range of values of a SV H and oifs shown in these windows encompasses most of the physically allowed 
range, namely from zero to J\f (in which case every incident atom induces a reaction). In both cases the region 
labeled saturated shrinks with increasing values of the efficiencies a sv n and aFs, and eventually disappears when 
a SV H + 2 aFs = 07at (see ec f @)- The reason is that both processes involve release of hydrogen from the grain 
and, consequently, when a SV H +2a Fs = O the grain can be sustained unsaturated irrespective of the thermal and 
quantum release rates. For example, the dotted and dash-dotted curves in windows (c) and (d) correspond to a case 
where the grain is unsaturated in the entire (— lnT^/t'o, E e /kT g ) plane. The panels show that the H 2 formation 
efficiency decreases with increasing a sv n and increases with increasing ap s - Because case (c) involves release of 
H atoms by direct collisions, the efficiency curves become independent of E e for values larger than that at which 
the collisional ejection rate exceeds the thermal evaporation rate, that is, Th < cx sv hS. Likewise, case (d) involves 
collisional formation followed by ejection of H 2 molecules and, therefore, the efficiency curves become independent of 
^ 2 h/^o for values such that T 2 H < aFsS. 


6 . H 2 from Surface Formation plus Pair Evaporation and Prompt Reaction 


We now generalize the above analysis to grains whose surfaces can form and dissociate H 2 while continuing to 
assume no H 2 formation/dissociation in the bulk. Consequently, the rate constants tp s and tjj) assume non¬ 
trivial values and the definition of cxfs is modified. The dimensionless coefficients 05 and ae acquire non-zero values 
but the previous analysis remains basically unchanged. In particular, the rate coefficients 03 and 04 and the expression 
for 9 sa t remain unaltered. H saturation in the bulk still occurs when <24 > gb irrespective of the values of the other 
parameters. The presence of H 2 formation on the surface does not change the characteristic release rate of H from 
the bulk. (This would not be true in the presence of H 2 formation in the bulk, because then hydrogen can be released 
from the bulk also in molecular form.) We proceed, as before, by setting 04 = 0. For q = 0, (no ejections of H atoms 
following collisional dissociation on the surface) eq. & which gives the condition for H not being saturated, and 
eq. (^||), which gives the efficiency, are modified as follows: 


T2 H r 2H + t F \ 


= r 2 # + 


2 

t dI + r ff 2 


(34) 

(35) 


Since the rates of surface formation and pair evaporation are each quadratic in 9h 1 the previous analysis is only 
altered by the change in the effective H 2 formation rate coefficient. The indicated change above clearly shows that 
the efficiency will increase when H 2 can form on the surface but the increase is qualitatively important only when pair 
evaporation and prompt reactions are ineffective by themselves {S{T 2 h / 9s^ 2 h + ^Fs/^h) << 1 ) while the additional 
process is sufficiently rapid ( Sip s /g s Tf I >1). The form of tpl depends on the ratio of surface dissociation rate to 
evaporation rate for H 2 . If evaporation is more rapid then tp —> tp^ and nearly every surface-formed H 2 escapes as 
soon as it is created. Efficient molecule formation requires St~f} g /g s T 2 H > 1. On the other hand, if dissociation is more 
rapid then H and H 2 obey a Saha-like equilibrium on the surface and the molecule escape rate is the evaporation rate 
of the H 2 fraction. Efficient formation requires STH 2 z 2Ed ^ kTg /9 sF 2 h < 11 - 

The formation of H 2 on the surface provides a source of molecular hydrogen on and in the grain and one must 
worry about the possibility of H 2 saturation in addition to H saturation. The region in the space spanned by the 
parameters ai through 07 in which the grain is unsaturated by any component is determined, in terms of the solution 
of eqs. (^5|)-(pg|), by the set of inequalities, 
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d H (ai,a 2 ) + Oh 2 («ii a 2 , « 5 , a 6 ) < 9s 

n ff (ai,a 2 ,a 3 ,a 4 ) + «h 2 (ai,a 2 ,a 5 ,a 6 , 07 ) < 9 b- (36) 

This also defines the regime of unsaturated grain in the space of physical parameters through the dependence of 01-07 
on them. 

Rather surprisingly, it is possible to give an approximate compact analytic treatment of the above conditions. We 
have the limiting solutions for 9h 


J 1 /ai for of >> 4a 2 
\ l/o ^ /2 for of << 4a 2 


and the condition that the grain (both surface and bulk) be unsaturated is 


9 2 h « min I —, 


a 5 \ 1 + °6 


9b — 0-4: f 9b ^4 


0507 \a 3 + 0607 


(37) 


(38) 


Terms involving g s ( gi,) are the surface (bulk) constraints. (As discussed previously, if g b < 04 the bulk must become 
saturated.) These inequalities guarantee that the regions of saturated and unsaturated grains will be identical in 
shape to those shown in window (a) of figure [i] in which the axes are relabeled, viz. 


E e ^ _ ln Sai 
kT g v 0 


r 2H 



(39) 


and the numerical values for the horizontal and vertical lines are determined from the above list. 
It is also straightforward to show that the limiting forms for efficiency for q = 0 are 


2a p 


£ = 


( rH+2op a s 

l-( 


1/2 


\2S(r 2 H+ipl) y 


for a\ » 4a 2 
for af << 4a 2 


(40) 


When <;^0, the efficiency is lowered by the amount 2qf3 s 9n 2 ■ 


7. H 2 Formation in the Bulk 


Next, we wish to examine the conditions under which effective H 2 formation may take place in the bulk, and explore 
the different regimes in parameter space in which grain saturation by H and/or H 2 occurs. Since the inclusion of 
quantum pair evaporation complicates the analysis considerably, and can only lead to a higher H 2 formation efficiency 
anyway, we ignore it in the present example. We also ignore H 2 formation and dissociation on the surface (but keep 
the prompt reaction process) in order to simplify the analysis further. The complementary case in which H 2 formation 
is controlled entirely by the surface has been studied in the previous example. 

Consider a grain for which H 2 formation in the bulk is allowed but H 2 formation and destruction on the surface as 
well as quantum pair evaporation are negligible. We set tp] = Sf s = Tp] = T 2 h = 0 in eqs. (20)-(p3|). One then 
finds, 


1 - bin H — b 2 n 2 H = 0, 


with 


^ _ ( 2 oiFb _j_ / gsTsbH_\ Th/S + 2 a Fs 

\1 +r] \gbnsH J l + T sbH (T H + ZotFsS) 

, 2 

O 2 = -- T 

g b tFbS( 1 + g) 

and we have introduced two auxiliary variables to simplify the notation: 


(41) 


(42) 
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(43) 


_ ffb(l + T sbH 2 T H 2 )T b sH 2 
g s TsbH 2 ^ H 2 TDb 

1 _ Sh_ __ 

r S 5(1 + T s bH (Th + 2aF S S)) 

It can be easily verified that in the absence of H 2 formation inside the grain, which corresponds to the limit tW = otFb — 
0 , the solution of eq. pl|), given by nn = l/&i, coincides with the solution of eq. ( p5| ) and (^ 8 |) for T 2 h = 0 , td s = 0 
as required. Note that the quadratic coefficient b 2 yf 0 only if H 2 can form by thermal processes within the bulk and 
its effect is to lower nn ■ The complexity of the linear coefficient &i arises from all the linear processes that alter the 
distribution of H (and H 2 ) in and on the grain including evaporation from the surface (T//, T// 2 ), transfer between 
bulk and surface ( T s bH, T~b S H, T s bH 2 , T~b S H 2 ), and collisional H 2 formation by particles incident on the bulk and on the 
surface (a.Fb, of s ). 

The quantity 1/r is the equilibrium fraction of the incident H atoms available to the bulk, r lies in the range 
between 1 , for Sh — S (as in the case of homogeneous injection for which Sh = g s S and Sh = gbS ~ S , provided 
that the grain is large enough), and 1 + TabH^H + 2 cxf s S) for Sh = 0 (which corresponds to a situation where the 
ambient gas is cold such that atoms impinging on the grain stick to the surface and diffuse into the bulk). In the latter 
case, nn approaches zero when the surface to bulk exchange time t s i,h becomes much longer than the evaporation 
and prompt reaction timescales since the only source of bulk atomic hydrogen is the flux from the surface. 

The coefficient g is the ratio of dissociation rate and net release rate of molecular hydrogen in the bulk. If rnb is 
much longer than the H 2 release time, viz., 77 << 1 , almost every newly formed molecule will escape from the bulk 
before it dissociates and, consequently, the loss rate of bulk atomic hydrogen will depend on the molecule formation 
rate (thermally activated and collisional) alone. When g » 1, only a fraction g~ l of the newly formed molecules 
leave the bulk before they dissociate, so the loss rate of H in the bulk is reduced by a factor g _1 . 

The surface occupation probability of H satisfies the equation, 


Oh = b 3 n H + b±, 


where eqs. (j||) and (|||) have been employed, and where 


b 3 


64 


_ g s (r s bH/U sh) _ 

9b[ 1 + TsbH^H + 2of s 5)] ’ 

_ T sbH<SH _ 

1 + TsbH^H + 2 CXFsS) 


(44) 


(45) 


The parameter 63 represents the net contribution from the bulk to surface by exchange, whereas 64 represents the 
contribution from direct injection of H atoms on the surface. When t s 5 / 4 ( 1 " H + 2 aFsS) » 1, the loss of hydrogen 
atoms from the surface is dominated by prompt reaction, evaporation and direct ejections rather than transfer to 
the bulk. Then (Th + 2 cxfsS)6h — gsT^nip-H / gb) + Sh ; that is, the loss of surface H atoms by evaporation, direct 
ejection, and prompt reaction is balanced by injection into surface sites plus transfer of H atoms from the bulk. In 
the opposite limit, ( On/gs ) — (jsbH/Tb S H)[nH/gb) + t^hSh / 3 s, which for sufficiently small surface injection rate 
(jsbHSn — 0 ) approaches the equilibrium distribution (see discussion following eq. Q). 

The bulk and surface concentrations of molecular hydrogen depend on three additional parameters: 


b 3 = 


TDb ( g 


g b t F b \ 1 + V 


be — aF b STF>b 

b 7 = 1 


V 


1 + g 


gT DbT h 2 

In terms of these parameters, the solution for nn 2 takes the simple form 

tih 2 = b 5 n 2 H + b e uH, 

and that for 9h 2 , 

0h 2 = b 7 riH 2 ■ 


(46) 


(47) 


(48) 
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The last equation reflects the fact that the only source of H2 molecules on the surface is the flux of H2 from the bulk. 
This is due to our neglect of H2 formation on the surface. 

In general, the regime of validity of the solution presented above is defined by: 

n H (b 1 ,b 2 ) +n H 2 {bi,b2,b 5 ,b 6 ) < g b , 

dn{bi 1 b2^b'i, 64) + 9 H 2 (bi,b 2 ,b 5 , be, by) < g s . ( 49 ) 


We now note that the substitution —> a*, the interchanges n <-> 6 and g b <-> g s yield exactly the same set of equations 
analyzed in the preceding section. Thus, we have the limiting solution ur 


( 1 /bi for b\ » 4 & 2 
UH \ l/b ^ 2 for b\ « 4 & 2 

and the condition that the grain (both surface and bulk) be unsaturated is 

2 


n~ H « min 



gs b<\ { g s 64 


b$by \be + beby 


( 50 ) 


( 51 ) 


When 64 > g s , the injection rate of hydrogen on the surface exceeds the rate at which H atoms are lost from the 
surface by evaporation, prompt reaction, and transfer to the bulk, as can be seen from eq. (phi), and, therefore, surface 
saturation will occur regardless of the value of hr. 

In contrast to the problem of H2 formation studied in the previous section, there is no simple relationship between 
b\ and b 2 and the fundamental microphysical parameters (e.g. E e /kT g and Gr )■ This is already apparent in the 
original expressions for b\ and b 2 which involve complex relationships between all the formation and release processes 
for H and H2. Nonetheless, the above inequalities faithfully convey the conditions for surface and bulk to remain 
unsaturated. 

We define the efficiency of H 2 formation as in the previous example: £ = 1 — TrOr/S. By employing eqs. (pi] ) and 
(pl|), we can express £ as 


£ = 1 — ThS 1 [637144(61, b 2 ) + 54 ] = 1 — T# S 1 



-bi + 




( 52 ) 


The efficiency is independent of the parameters 65 through by, as it measures only the ratio of escape rate (thermal 
evaporation plus direct ejections) to interception rate by the grain of H nuclei. It does not tell us anything about the 
values of bulk and surface concentrations of molecular hydrogen. The latter are determined by the rate coefficients 
65, be and by, which involve ratios of H2 formation rates to dissociation, release, and bulk to surface exchange rates, 
as seen from eq. (p6|). 

To elucidate the solution further we consider the limit a Fs = 0 and Sh = 0 (r = 1 ). The latter approximation 
implies 64 = 0 . We define 


x F = l/t Fb S, 

_ g s (T s bH/T b sH)(rH/S) 

<7i>(l + TsbH^H) 


( 53 ) 


x F is the thermally activated H2 formation rate in the bulk in units of the injection rate per site, S. xr is the 
net release rate per bulk H atom of atomic hydrogen from the bulk measured in units of S. For h » 1 ) 

xr — (g s /g b ){ T bsHS )^ 1 , in which case H release from the bulk is limited by the rate at which H is transferred to the 
surface. For r sb H^ h « 1 , xr — {g s /gb){T sb H/x b sH)(rH / S) and the net release of H is limited by evaporation from 
the grain (g s T sb H / gbX bs H is the equilibrium ratio of surface and bulk H concentrations). 

Eq. © may be recast as a relation between the formation rate and the release rate for fixed hr 


xr = 1 rijj 1 


2 a Fb 


)- 


2 ur 


-x f . 


Likewise eq. 


1 + vJ 9biX + v) 

may be used to find contours of fixed efficiency in the (xr,x f ) plane 


( 54 ) 


g b [£{ 1 + y)x 2 R - 2(1 - £)a Fb x R ] 
2(1 — £) 2 


( 55 ) 
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Finally, once 64 , 65 , be, and 67 are fixed, the critical value of tir for which the grain becomes saturated in or on can 
be computed from eq. (|49|). 

Contours of hr in the (xr,xf) plane are exhibited graphically in figure || (values of hr label the curves), for 
&Fb = 0, 77 = 0 (solid lines) and g = 10 2 (dotted lines). When the release rate is large, then hr is small (near the 
top). When the rate is small and the formation rate is small, then ur is large (bottom right); while if the formation 
rate is large, then it is small (bottom left). Another way to view the results is as follows. The vertical portion of the 
curve (independent of xr) corresponds to the regime where uh is limited by H 2 formation; the horizontal portion 
(independent of xf) is the regime where hr is limited by release from the bulk. As g increases the curves move 
leftward. The reason, as already explained above, is that for g > 1 the rate at which hydrogen leaves the bulk in 
molecular form diminishes by g _1 and, therefore, shorter formation time (i.e., larger xf) is required to maintain tlr 
at the same value, as also clearly seen from eq. ([i4j). Also shown schematically in figure is the curve separating 
the regimes of saturated and unsaturated grain (dashed line). The location of this curve depends on the values of the 
remaining parameters (fo - 67 ) through eq. (|49|). The corresponding critical value of ur can range between zero for 
64 > g s and ~ gb- 

The quantitative range of uh is easily derived: For a given injection rate, uh cannot exceed n 0 = (1 + g)/2aFb , 
when capture of H atoms in the bulk is balanced by release of H 2 molecules produced by direct collisions. For this 
value eq. (|54|) yields xr = 0, Xf = 0. Sustaining hr at values smaller than n 0 requires, in addition to collisional 
formation of H 2 molecules, either thermally activated H 2 formation (xf > 0), or H release (xr > 0). For uh « n 0 
the latter processes dominate over collisional H 2 formation which ultimately becomes negligible. 

Labeled efficiency contours are presented in figure ||for g = 0, aFb = 0 (solid lines) and ctFb = 10 4 (dashed lines). 
As seen, for large Xr (large release rate) the efficiency curves are power laws with slope 2, reflecting the dominance 
of thermally activated over collisional H 2 formation. As Xr approaches 2(1 — £)apb/£ , the desired efficiency can be 
achieved by collisional H 2 formation alone and, hence, Xf —> 0 . 


8 . Saturated Grain 

What happens when the grain becomes saturated by either atoms or molecules? Presumably, the saturation of 
the grain affects the embedding energies and diffusion times considerably. Unfortunately, the dependence of those 
parameters on filling factors is poorly understood. Nevertheless, we can still draw some general conclusions. Consider 
a clean grain having a temperature T g , initially devoid of H atoms and H 2 molecules, embedded in a gas of hydrogen 
atoms. Suppose first that molecular hydrogen formation in and on the grain is energetically unfavorable and/or 
inefficient (i.e., tpi —> 00) and that the mean evaporation time of an H atom, the prompt-reaction time and the mean 
surface recombination time are all much longer than ti n . Then after a time ~ A fti n the grain will become saturated by 
H atoms. The evolution of the system from this point on will depend on the modified parameters. If the embedding 
energy of H atoms is not altered significantly, then a newly incident atom will either eject bound H atoms from the 
grain as a result of collisions provided that its energy is sufficiently large, push bulk atoms to the surface, or will 
stop inside the grain and then diffuse to the surface and stick. If the yield for the first process is much smaller than 
unity, then layers of atomic hydrogen should form on top of the original grain surface. Our expectation is that this 
will lead to a dramatic increase of the surface recombination rate and, consequently, effective H 2 formation, since the 
atom-surface binding energy declines sharply with distance from the original grain surface. 

Next, consider the possibility that molecular hydrogen formation in or on the grain is effective, and that the mean 
evaporation time of a molecule is much longer than 2t,„. The grain will become saturated by H 2 molecules after a 
time ~ 2J\fti n . A newly injected atom will then either break up a bound H 2 , ejecting one of the H atoms from the 
grain in the process, or be pushed from the bulk to the surface and then evaporate or stick, depending on the modified 
surface binding energy. If it sticks to the surface for a time longer than about U n , it will eventually recombine with 
a second injected atom to form a molecule (under astrophysical conditions, the surface migration time is anticipated 
to be much shorter than ti n ). This will lead to the growth of layers of molecular hydrogen until the evaporation rate 
of atoms or molecules exceeds the value required for steady-state. 

Alternatively, it could be that when the filling factor becomes sufficiently large, a few per cent say, the embedding 
energies of H and H 2 are sufficiently reduced to allow the evaporation rate of either atoms or molecules to balance 
the flux of incident particles. The efficiency of molecular formation in this scenario would depend on the value of the 
modified parameters at steady-state. Our analysis should be applicable provided that the hydrogen has not undergone 
a phase transition to some ordered phase at these densities. 
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V. APPLICATION: H 2 FORMATION IN CARBON 


As an example we consider H 2 formation in and on carbon grains. For simplicity, we take the distribution of H 
atoms in the gas phase to be isotropic and monoenergetic with respect to the grain, and denote by n(H) and e the 
number density and energy (measured in units of eV) of the H atoms. The velocity of an H atom relative to the grain 
is then v = 1.4 x 10 6 (e/eV) 1 / 2 cm s -1 , and the flux impinging on the grain is n(H)v ~ 4 x 10 5 (e/eV) 1,/2 n(H) cm -2 
s -1 . In terms of the transmission coefficient T(e), defined in the text preceding eq. ([l]) , and the characteristic grain 
size a = a_ 5 l 0 -5 cm, we can approximate the injection rate as 

<^3x 10- 4 n(H)(e/eV) 1 / 2 ol 5 T(e) s" 1 , (56) 

where the grain’s geometrical cross section has been taken to be of order 47ra 2 (ignoring possible focusing effects). 
Let d-8 be the average interatomic separation of C atoms in A. Then the total number of sites is roughly A f ~ 
4 x 10 9 (a_ 5 /d_ 8 ) 3 and the relative number of surface sites is of order g s ~2x 10~ 3 (a_ 5 /d_ 8 ) _1 . Thus, 

S = (A/fin) -1 — 7.5 x 10” 14 n(H)(e/eV) 1/2 aIgd'L 8 T(e) s _1 . (57) 


Hydrogen adsorption on carbon takes place in different types of sites with different activati on .energies. The 
experimental values of the activation energies for adsor ption lie between about 0.3 and 2.2 eVuO. The heat of 
adsorption (or equivalently the embedding energy E e , see § III D and figure^) has not been determined experimentally 
yet. However, .-several attempts have been made to calculate it for the hydrogen graphite system using quantum 
chemical codesE§E3. The results of such calculations are highly uncertain, but seem to suggest the following: i) the 
most favorable site for chemisorption is the one which corresponds to the H atom being directly above a C atom. 
The embedding energy lies in the range between ~ 0.3 and 1.14 eV, implying E c between 1. 94 an d 1.1 eV. ii) the 
activation energy for H pair recombination (E r in figure [)]) is about 0.4 eV. As mentioned in § [II,L| , the detection of 
highly excited H 2 molecules formed on carbon surfaces at low temperatures led Schpwnann et ahEJ to conclude that 
the embedding energy should not exceed ~ 0.3 eV, consistent with the results of ref.ca. For E e in the range 0.3 — 1.1 
eV and E r = 0.4 eV, we obtain an activation energy of chemisorption in the range between 2.34 (note that this implies 
that there should be an activation barrier of ~ 0.1 eV for adsorption of atomic hydrogen) and 1.5 eV, respectively, 
consistent with the experimental values mentioned above. 

Taking E r = 0.4 eV, eq. © yields for the penetration factor at zero temperature Gh — 30 (z 0 /A), and eq. 
gives T t — 150 (z 0 /A) K for the tunneling temperature, where z a is again the width of the barrier. In the absence 
of collisional processes eq. ([n]) with 9 sat = g s implies that the surface will be maintained unsaturated at T g = 0 
provided that 


ln[n(H)(e/eV) 1/2 d?_ 8 T(e)] < 54-30 (z 0 /A). 


(58) 


As indicated above, we are ignoring the possibility that bulk saturation alters fundamental grain parameters and 
assume this condition to suffice even if the bulk is saturated. (We show that bulk saturation is likely below.) Adopting 
z 0 = L4, d = 3 A, T(e) — 1 for illustration, we find that quantum pair evaporation will keep the surface unsaturated 
as long as n( H)(e/eV) 1//2 < 2.5 x 10 9 cm -3 . Increasing the grain temperature will result of course in a smaller surface 
concentration. Note that the condition (|5|) is independent of the grain size since both the incoming flux of H atoms 
and the outgoing flux of H pairs are proportional to the surface area. 

If molecules form only by pair evaporation from the surface, the efficiency is given by eq. (^3|). Setting the collisional 
rates to zero and taking £ = 1/2 in eq. (|3^), we find that for the same choice of z a , d and T(e), the efficiency will 
exceed 1/2 if 


T a < T crt = 


10 i (E e /leV) 


40 — ln(n(H)(e/eV) 1 / 2 )/2 


K. 


(59) 


For example, for the value of E e inferred by Schermann et al., viz., E e = 0.3 eV, the last equation gives T crt between 
about 75 and 100 K for n(H)(e/eV) 1 / 2 between 1 and 10 7 cm 3 , respectively. It is important to note that the efficiency 
£ defined in eq. (( 2 J) is essentially the probability that an H atom intercepted by the grain will be released as part 
of a molecule. The overall H 2 formation efficiency is the product of £ and the transmission (sticking or penetrating) 
probability of H. 

If collisional displacement of H from bulk to surface does not occur and if molecule formation in the bulk does not 
occur then we find that the bulk will be saturated by H in steady state. Our reasoning follows. The solubility and 
diffusivity-iof hydrogen in the bulk are uncertain. The reported values of the diffusion energy range between 0.5 and 
4.3 eVE3’E3. This is most likely due to the broad spectrum of binding sites. The deepest sites are the trapping sites 
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thought rta-be associated with unsaturated carbon atoms on the edges of microcrystals constituting the amorphous 
materialOCJ. The trap energy is of order 4.3 eV, comparable with the bond energy between H and C in typical 
hydrocarbons. The mobility of hydrogen at low concentrations is presumably controlled by the deep sites. In fact, 
implantation experiments appear to suggest that atomic hydrogen may be highly mobile in the absence of lattice 
damages^. Astrophysical grains are likely to be damaged. However, once the deep sites are filled the mobility of 
additional atomic hydrogeruin the bulk could be large. The best value for the solution energy, according to ref.ua, is 
that given by Atsumi et alO: E s = —0.2 eV (cf. figure |l| case iv), and for the diffusion energy is Ed = 2.8 e V. The 
energy difference between the surface and bulk ground states is given by E e — E[H] b = E c — E s (cf. § jlll D ). With 
the values of E c discussed above, E s = —0.2 eV implies that this energy gap lies in the range between 1.3 and 2.1 eV. 
The condition on the maximum equilibrium abundance to prevent saturating the bulk is very strict on account of the 
large energy gap: 0 sat = g s e~^ Ec ^ Es ^ kT s (cf. eq. (30| ). If molecule formation on the surface is efficient we require 


ln[n(H)(e/eV) 1/2 cP 8 T(e)] < 54 - 30 (z 0 /A) min 1^ - 


3 x 10 4 (E c -E, 


T n 


1.3eV 


or if H evaporation dominates 


\n[n(R)(e/eV) 1/2 d 2 _ 8 T(e)} <54- 


2.8 x 10 4 / E e + E r — E< 


T n 


2.44eV 


(60) 


(61) 


Neither condition can be satisfied for typical astrophysical conditions so H saturation is favored. 

It is important to recognize that the surface to bulk exchange rate has a large typical activation energy A E = 
Ed + E s — E c . For Ed = 2.8 eV, we obtain A E between 1.5 and 0.7 eV. Given these estimates we anticipate the 
surface to bulk exchange rate to be much smaller than H or quantum pair evaporation rate. Bulk saturation is 
achieved, however, if high energy particles penetrate the lattice. 

This relatively simple picture for graphite could be complicated by the following additional chemical pathways. The 
bulk may avoid saturation in steady state if sufficiently rapid H 2 formation can occur in the bulk followed by rapid 
bulk to surface exchange of H 2 , or in the presence of rapid collisional transfer of bulk H to the surface. 

Whether H 2 formation can take place in the bulk is unclear. It can be argued, based on the foregoing discussion, 
that if H 2 is stable in the bulk then the activation barrier for bulk to surface exchange of H 2 molecules is large (it must 
exceed that for H atoms; c.f figure ||) so that H 2 saturation is inevitable for the same reasons that H saturation is. In 
both cases, saturation of the deep binding sites could alter the effective interaction potential. By simply occupying 
the deepest sites, hydrogen within the grain may force newly added atoms or molecules to move amongst more weakly 
bound sites. Alternatively, completely new hydrogen interaction potentials may be relevant at high occupancies. In 
any case, the high abundance of H and/or H 2 can only be altered from that predicted previously if rapid migration 
to the surface of H or newly formed H 2 becomes possible. In fact, in refsEj it has been argued that there is strong 
indication of H 2 formation and migration in some forms of graphite and amorphous carbon films saturated with 
hydrogen. There is insufficient data to make a quantitative estimate of the rate for the latter process. It is likely, 
though, that once the bulk becomes saturated, then a newly intercepted atom in the bulk will either diffuse quickly 
to the surface where it evaporates or undergoes quantum pair recombination, as discussed above, or reside in the 
bulk long enough to recombine with another incident atom to form an H 2 molecule that will subsequently be ejected 
from the grain. The rates for these processes depend on the barriers of the shallow sites in the saturated bulk and 
the modified activation energy for bulk to surface transition. In all these more complicated scenarios, H 2 formation 
increases compared to the simple estimates based solely on pair evaporation. Thus, we conclude that at low enough 
grain temperatures effective H 2 formation should take place in or on these putative carbon grains. 


VI. SUMMARY AND CONCLUSIONS 

In this paper we have examined the conditions required for effective molecular hydrogen formation in and on solids, 
with particular emphasis on astrophysical grains. 

Much of the earlier work has been concerned with the catalysis of H 2 formation by surfaces of grains in molecular 
clouds, where the grain and gas temperatures are very low. The key processes determining the rate of H 2 formation 
are then: sticking and retention of the gas-phase atoms, diffusion of H along the surface, recombination with another 
H atom, and ejection of the newly formed H 2 molecule. In these scenarios the chemisorption sites are taken to be 
saturated and hydrogen atoms^and molecules are weakly bound to physisorption sites with binding energies on the 
order of a few tipes 10 -2 eV0. This path may apply to grains coated by overlayers of molecular material such 
as amorphous iceua. The evaporation time of H from these weak binding sites is very short even at the low grain 
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temperatures envisioned- but deeper (impurity) sites may retain a fraction of the H atoms for long enough time, 
and act as reaction sitesB. The migration time of, a newly adsorbed H atom to a reaction site is typically very short 
compared with the residence timeE3 (but c.f. ref.Ea), so that the factors limiting the H 2 formation rate are the sticking 
probability and retention of H. 

The mechanisms controlling the rate of molecular hydrogen formation on grains could be markedly different un¬ 
der harsher conditions, e.g., in the vicinity of an astrophysical shock. Firstly, the grain may be exposed to fluxes 
of energetic particles or intense UV/X radiation, and this may lead to surface cleaning, surface reshaping, lattice 
modification and so forth that create chemisorption sites. Secondly, the gas-phase atoms may be energetic enough to 
penetrate deep inside the grain and, under certain conditions, to drive direct collisional reactions. H 2 formation may 
then occur, (1) in the bulk via thermal activation, (2) in the bulk by direct recombination with incident H atoms, 
(3) on the surface by thermally activated recombination of chemisorbed H atoms, (4) by direct recombination of an 
incoming atom with a chemisorbed atom, and (5) through pair evaporation. In case 4 a fraction of the newly formed 
H 2 molecules may be ejected immediately from the grain, a possibility discussed earlier in ref.c3, and referred to .hepe, 
following him, as prompt reaction. Chemisorption of hydrogen on carbon has been studied by several authorsr 18 H 49 l, 
who calculated interaction potentials and chemisorption energies of atomic hydrogen on ideal graphite surfaces using 
quantum chemical techniques. There remains some variation in these results, so we have treated the value of the 
binding energy as a parameter. 

The efficiency of H 2 formation depends, quite generally, on the characteristics of the solid including the following 
energy scales: the ground state energies of H and H 2 in the bulk, surface and vacuum; the energy barriers for site 
to site diffusion, for surface to bulk exchange, for H 2 formation in the bulk, on the surface and for pair formation 
above the surface. The efficiency also depends on the grain temperature and the flux and energy (particularly if high 
enough to drive collisional reactions) of H nuclei impinging on the grain. It is important to remember that the energy 
levels and activation barriers themselves depend on the density of hydrogen in and on the solid by virtue of the self 
interaction of hydrogen atoms and molecules. These values could be altered substantially when the grain reaches 
saturation. 

In this work we have explored various recombination reaction pathways that may operate under a wide range of 
conditions. By constructing simple analytic models, we identify the regimes in the space of physical parameters that 
correspond to effective H 2 formation and grain saturation by H atoms and H 2 molecules. Although our model is 
general we have focused on the case of “clean” grains with only one type of surface site and one type of bulk site. So, 
our examples do not apply to the classic molecule formation mechanism explored by Hollenbach and SalpeterEil which 
envisioned two surface sites with different binding energies. We also assume that hydrogen diffusion in the bulk and 
on the surface is rapid enough to render the bulk and surface occupation probabilities homogeneous. The processes 
incorporated in our model include: thermally activated and collisionally induced exchange of H and H 2 between the 
bulk and the surface, H 2 formation in the bulk and on the surface by thermal activation and direct recombination 
with incoming atoms, pair evaporation, prompt reaction, thermal evaporation and collisional ejections and a variety 
of dissociation reactions in and on the solid (by thermal activation, by quantum tunneling and by direct collisions). 

Despite the large number of physical parameters involved, the solutions for the H and H 2 bulk and surface con¬ 
centrations and the H 2 formation efficiency are found to be highly degenerate, and depend only on a relatively small 
number of dimensionless parameters. In the following, we briefly summarize the main results for the three distinct 
cases we have considered: when H 2 forms from H atoms leaving the surface, when it also forms on the surface, and 
finally when it forms in the bulk. 

If there are no bound H 2 states in the bulk and on the surface, H 2 formation occurs via pair evaporation or prompt 
reaction. Two cases are of interest, (a) In materials which have positive chemisorption energy, as appears to be the case 
for carbon, pair evaporation can proceed quantum mechanically even at zero temperature. Our result is that £ > 0.5 
when S{{T 2 H/gs^ 2 H)+^ a Fs/^ h) > 1- In the absence of collisions this is simply — ln(T 2 H/v o ) — 2E e /kTg—\n(S/g s v 0 ) < 
0 so that effective H 2 formation is “guaranteed” as long as the loss rate of atomic hydrogen by evaporation is less than 
the quantum pair evaporation rate. This is essentially a condition on the maximum grain temperature. (The threshold 
grain temperature does depend on the rate at which H nuclei are intercepted by the grain but the dependence is 
logarithmic and therefore weak.) If the grain temperature is large, the surface density of H is small and the efficiency 
of formation of H 2 , which is a process quadratic in the surface density, is small, (b) For materials with negative 
chemisorption energy, binding to the surface of H is energetically favorable so that surface densities would saturate 
at typical grain temperatures. Pair evaporation and H evaporation both require thermal activation but generally 
molecule formation will proceed with more rapidity simply because H 2 is bound with respect to H in the vacuum. 
(We are assuming that the barrier to H 2 formation is less than the embedding energy of H.) For these materials, 
the efficiency of molecule formation will be high. However, since the overall loss rate is exponentially sensitive to 
temperature these materials will become saturated if the grain temperature is too low. Saturation of the surface 
occurs when 2T 2 h + T h < gf 1 S- Whenever grain saturation of the surface occurs it is possible that the important 
parameters describing the grain are modified. 
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For both materials above, collisional ejections of H atoms and prompt reaction may alter the results; the former 
releases atomic hydrogen from the grain and reduces the H 2 formation efficiency, whereas the latter leads to effectively 
higher surface recombination rate and, therefore, enhanced efficiency. The rates of both processes are independent 
of the grain temperature (to the extent that aFs and a sv themselves do not depend on T gr ). As figure illustrates, 
these processes can significantly alter the regimes of effective and ineffective H 2 formation and also grain saturation. 

When H 2 can associate/dissociate and adhere to the grain surface (but not the interior) the above situation is 
more complicated. H 2 formation on the surface always increases the formation rate by pair evaporation alone [this 
is mathematically made clear by the simple modification needed for eqn. (|3l|)1. The efficiency of molecule formation 
is increased. This qualitatively modifies our previous example only when formation would be otherwise inefficient. 
Efficient formation requires tp] » g s T 2 H /S, where tp] is the effective formation rate defined explicitly in eq. (|2£|). 
The limiting cases were discussed and correspond to (1) H 2 formation followed by immediate evaporation from the 
surface and (2) evaporation of the surface equilibrium H 2 fraction. Numerically, efficiency contours are given by the 
modified form of eq. (|ts|) . 

In so far as saturation is concerned, it is more difficult for H to saturate the surface. At the same time, surface 
formation also introduces the possibility that H 2 may become saturated on the surface if the release rate from the 
grain is not fast enough. The detailed conditions are given in eqn. (|38|). 

H 2 production in the bulk requires the existence of bound H 2 states inside the solid, as depicted in figure [|. 
It is unknown which substances allow such bound states. If bound states exist, the process has many similarities 
to H 2 formation on the surface. Likewise, it increases the net efficiency. We have calculated the effect of bulk 
formation independent of the surface and pair evaporation processes. The H and H 2 occupation probabilities as 
well as the formation efficiency depend on the net rate of release of atomic and molecular hydrogen from the bulk 
and, consequently, involve the exchange rates between bulk and surface. The full solution for H and H 2 occupancies 
and the efficiency depends on 7 dimensionless parameters but the bulk concentration depends on only 2 . Additional 
simplifications are possible if the transfer time of H from surface to bulk is very short or the injection rate of H into 
bulk sites exceeds that into surface sites. The latter case is particularly relevant for homogeneous injection, the limit 
achieved when high energy particles penetrate the lattice. We show that the bulk concentration of H may be simply 
expressed as a function of xr and Xf (figure ^]), where xp is the ratio of bulk H 2 formation rate to injection rate 
per site, and xr is the ratio of net release rate of H from the bulk to injection rate per site. The sensitivity of these 
results to the dissociation rate is also discussed (the ratio of H 2 dissociation rate and release rate of H 2 from the bulk 
is 77 ). The domains of H and H 2 saturation are described. 

H 2 formation within the bulk can proceed as a thermally activated reaction or by collisional impact of a high energy 
particle. If collisional H 2 formation is less important than thermally activated formation, then equilibrium is achieved 
between release and thermal formation. The efficiency depends on xp/x 2 R (see figure 0 and segment with slope 2). 
In the alternative limit, collisional formation achieves equilibrium with release (figure y with segment with slope 0). 
The location of the efficiency curves also depends on 77 . The efficiency is essentially independent of 77 when 77 << 1 
and decreases as 77 increases for 77 > 1 . 

As an application we have considered H 2 formation on carbon. This system is complicated by virtue of the broad 
spectrum of binding sites exists in the solid. Despite this complication we were able to draw some general conclusions 
using results of simulations and data available in the literature. In particular, 1) efficient molecule formation requires 
grain temperatures smaller than about 100 K, 2) pair evaporation will keep the surface unsaturated even at very low 
grain temperatures if 77 (H) (e/eV) 1 / 2 < 2.5 x 10 9 cm -3 where ?z(H) is the number density of gas phase atoms and e 
is their average energy, and 3) the bulk will be saturated in steady state unless sufficiently rapid molecule formation 
will take place in the bulk followed by rapid transfer of H 2 to the surface or if the incident atoms are energetic enough 
to eject bulk H to the surface at a large enough rate. 
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FIGURE CAPTIONS 


FIG. 1. Schematic energy diagrams of hydrogen in solids for which no bound H 2 states exist in the bulk and on the surface. 
The solid curve represents the interaction potential of a single H atom with the solid. The two dashed lines represent the 
interaction potential of an atom pair; one dashed line corresponds to a pair at fixed interatomic distance of order that of a free 
H 2 molecule, and the other to a minimum energy configuration. The chemisorption energy, E c , solution energy, E s , activation 
energy of chemisorption, E a , activation energy of pair evaporation, E r , and dissociation energy of H 2 in vacuum per H atom are 
indicated. The sign of E c and E s is defined by the direction of the corresponding arrows; positive if the arrow points upwards 
and negative if it points downwards. The four cases shown correspond to different possible ordering of the energy levels (see 
text for detailed discussions). 


FIG. 2. The same as figure |j], but in the case in which H 2 is stable in and on the solid. The thick arrows indicate different 
pathways for H 2 formation and release: [1] H 2 formation in the bulk, [2] H 2 formation on the surface, [3] pair evaporation, and 
[4] evaporation of H 2 from the surface. Each reaction may involve an activation barrier as shown schematically. 


FIG. 3. Domains of grain saturation, effective and ineffective H 2 formation for materials for which there are no bound H 2 
states on the surface and in the bulk, in the absence of (windows a and b) and in the presence of (windows c, d) collisional 
reactions. The independent variables are the log of pair evaporation rate per attempt frequency, T 2 h/vo, and the embedding 
energy in units of the grain temperature E e /kT g . 


FIG. 4. Conditions for effective and ineffective H 2 formation for specific materials. The dashed line corresponds to embedding 
energy E e = 2.55 eV, chemisorption energy E c = —0.3 eV, activation energy of chemisorption E a = 0.2 eV, and barrier 
penetration depth Gh > 7, and represents a typical transition metal. The dotted-dashed line corresponds to E c > 0, E e = 0.3 
eV, activation energy of pair evaporation E r = 0.4 eV, and Gh = 30, and represents carbon (see example in §|v|). 


FIG. 5. tih contours plotted in the plane defined by xf - the ratio of thermally activated H 2 formation rate in the bulk and 
injection rate per site, and xr - the net release rate per H atom of bulk atomic hydrogen in units of the injection rate per site, 
for ctFb = 0, rj = 0 (solid lines) and r) = 10 2 (dotted lines). Values of uh label the curves. The curve separating the regimes of 
saturated and unsaturated grain is shown schematically (dashed line). 


FIG. 6 . Contours of H 2 formation efficiency, £, for 77 = 0, ctFb = 0 (solid lines) and qf 6 = 10 2 (dashed lines). The efficiency 
curves are broken power laws (see eq. [^BJ). The change in slope reflects the transition from the regime where thermally 
activated H 2 formation dominates to the regime where collisional recombination of a bulk H atom with an incoming atom 
dominates. 
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APPENDIX: LIST OF SYMBOLS 

T g - Grain temperature 
T t - Tunneling temperature 

rix - Fraction of all sites occupied by species X in the bulk; (X =H, H 2 ) 

Ox - Fraction of all sites occupied by species X on the surface; (X =H, H 2 ) 

Af - Total number of sites 

gb - Ratio of the number of bulk sites and Af 

g s - Ratio of the number of surface sites and Af 

S - Total injection rate per site 

Sh - Injection rate per site into bulk sites 

Sh - Injection rate per site into surface sites 

£ - H 2 formation efficiency 

E[X } 1 - Ground state energy of species X; (X =H, H 2 , i = b,s, v ) 

E c - Chemisorption energy 

E e - Embedding energy 

E s - Solution energy 

E l d - Dissociation energy; (i = b , s, v) 

E r - Activation energy of pair evaporation 
E a - Activation energy of chemisorption 
Gh - Penetration factor for quantum pair evaporation 
ti n - Mean time for particle interception by the grain 
tpi - Thermally activated H 2 formation time; (i = b , s) 
toi - Thermally activated H 2 dissociation time; (* = b,s) 

tFsb ( tFbs ) - Time of H 2 formation in the bulk (on the surface) by thermally activated recombination of surface and 
subsurface atoms 

tDbs ( tDsb) - Time of H 2 dissociation in the bulk (on the surface) to produce a surface and subsurface H 

t s bx ( hsx ) - Time of thermally activated diffusion of a species X from the surface (bulk) to the bulk (surface) 

a.Fi - Probability per site of direct collisional recombination with an incoming H atom; (i = b , s) 

a sv x - Probability per site of collisional ejection of a species X 

Pi - Probability per site of collisional dissociation of an H 2 molecule; (i = b, s ) 

jx - Thermal evaporation rate of a species X from the surface 

Tjc - Total loss rate of a species X from the surface (sum of thermal evaporation and collisional ejection) 
r 2 h ~ Pair evaporation rate 

q - Fraction of collisional dissociations on the surface that results in the immediate loss of the H atoms 
Sfs - Fraction of H 2 molecules formed by collisional recombination on the surface that remain on the surface 
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